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Abstract 



^ ' The regular finite initial value problem at infinity is used to obtain regularity condi- 

O^i tions on the freely specifiable parts of initial data for the vacuum Einstein equations with 

non-vanishing second fundamental form. These conditions ensure that the solutions of the 

{3JQ| propagation equations implied by the conformal Einstein equations at the cylinder at spatial 

infinity extend smoothly (and in fact analytically) through the critical sets where spatial 

infinity touches null infinity. In order to ease the analysis the conformal metric is assumed to 

be analytic, although the results presented here could be generalised to a setting where the 

conformal metric is only smooth. The analysis given here is a generalisation of the analysis 

(■~^ ' on the regular finite initial value problem first carried out by Friedrich, for initial data sets 

l^-s , with non-vanishing second fundamental form. 

(N 

^ ! 1 Introduction 

oo . 

f ^ , The regular finite initial value problem near spatial infinity introduced in [11 provides an excep- 

tional tool for analysing the properties of the gravitational field in the regions of spacetime "close" 
to both spatial and null infinity. This initial value problem makes use of the so-called extended 

j^ ■ conformal Einstein field equations and general properties of conformal structures, and it is such 

JH I that both the equations and the data are regular at the conformal boundary. The formalism of 

the regular finite initial value problem has been used to analyse the behaviour of the conformal 
field equations on the critical sets where null infinity "touches" spatial infinity. The analysis in 
[TT] has shown that the solutions to conformal field equations develop a certain type of logarith- 
mic singularities at the critical sets. These singularities form an intrinsic part of the conformal 
structure. It could be the case that these logarithmic singularities do not affect the regularity of 
the rest of the spacetime — and in particular the smoothness of null infinity — but the hyperbolic 
nature of the propagation equations suggests otherwise. In 19 it has been shown that there is 
another class of logarithmic singularities forming at the critical sets. Subsequent generalisations 
of these last calculations carried out in [TSl [501 HI] seem to point out that stationary solutions 
play a prominent role in the discussion of the structure of spatial infinity. 

One of the crucial features of the formalism of the regular finite initial value problem at spatial 
infinity is that it allows to relate the behaviour of solutions to the conformal field equations at 
the critical sets with properties of the initial data. In particular, under the assumption of a time 
symmetric initial data set with an analytic conformal metric, it was possible to find conditions 
on the free initial data near spatial infinity which ensure that for solutions developing from from 
these data, singularities cannot occur. 

Given the need to systematise and gain further insight into the results for initial data with 
a non- vanishing second fundamental form given in [20[ I21j , in this article a generalisation of the 
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discussion of [TT] is carried out: it is studied how to construct conditions on free data of a class of 
non-time symmetric initial data so that their development does not have the type of logarithmic 
singularities discussed in [TT]. The understanding of the logarithmic singularities in P^fTSll^l^ 
is an outstanding open problem. 

In order to state the main result of this article, it is recalled that if one uses the so-called con- 
formal Ansatz to construct asymptotically Euclidean solutions to the vacuum Einstein constraint 
equations under the assumption that the initial data is maximal, then the freely specifiable data 
on, say, a compact mainfold S is given in terms of a pair of symmetric tensors {hij, $ij), where 
hij is the conformal metric of the initial hypersurface S and ^ij is /i-tracefree — throughout this 
article i, j, fc, . . . will denote spatial tensorial indices taking the values 0, 1, 2. Following the discus- 
sion in |5; , the tensor $ij is determined in a neighbourhood Ba G S oi the point at infinity i €z S 
by providing a real number A {the expansion), constant 3-vectors P* (the linear momentum), J* 
{the angular momentum), Q^ {the conformal momentum), and a complex function A which is C°° 
in Ba \ {i}- It can be seen that the real part of A gives a contribution, $y [A^-'^^], with properties 
similar to the contributions by A, P' and Q', hence it can be thought of as containing the higher 
mass multipoles of ^ij. On the other hand, the imaginary part of A gives a contribution ^ij [A^^] 
with properties analogous to J*. Accordingly, it can be interpreted as giving the higher angular 
momentum multipoles of ^ij . 

In terms of /i-normal coordinates a;' based on i, some technical assumptions will be made on the 
freely specifiable data. The precise details of these assumptions will be discussed in the main text. 
The conformal metric hij which in normal coordinates centred on i takes the form hij = —5ij +hij 
is taken to be in C'^{Ba) (the set of analytic functions on Ba), where in addition the technical 
condition 5"^^ hij = will be assumed. Furthermore, P* = and rA^^^rA^^^ G E°°{Ba), where a 
function / is said to be in E°°{Ba) if it can be written as / = f(^)+rf^^^ with Z^^^, J^^) g E°°{Ba), 
and r^ — (x^)'^ + (a;^)^ + {x^)'^ . This last assumption is a necessary condition for the solution 
of the Licncrowicz equation to be expanded in terms of powers of r only. In addition, it will 
be assumed that both A and hij are related in such a way that the solutions to the momentum 
constraint are of the form r^^^ij with $ij G E°°{Ba)- 

In order to express the regularity conditions product of the present analysis, a couple of further 
tensors will be required: 



hj 



~ j 



Duru - huDk^^'> 



4 



[A(^)]^6,;'Z?|,XZ|.)[A(^)], 



with Xiji-^^^''] the part of the second fundamental form constructed out of <i>,y[A*^''^-'] — obtained 
from *I'jj[A^^''] by multiplying by a conformal factor — , and r^, ^^'r ,Dk, respectively the Ricci 
tensor, the Ricci scalar and the Levi-Civita connection of hij, and eijk denotes the volume form 
on S. The tensor bij is the Hodge dual of the Coton(- York-Bach) tensor of hij. 

Theorem 1. For the class of data under consideration, the solution to the regular finite initial 
value problem at spatial infinity is smooth through the critical sets only if 

<^(A, •••Ai6jfe)(J) = p = 0,1,2,... 

<^(A, • • • AiC,A,.[A(«)]e^^^;)W = q = 0,1,2,.. . 

is satisfied by the initial data. If the above conditions are violated at some order p or q, then the 
solution will develop logarithmic singularities at the critical sets. 

In the above expression "^ denotes the operation of taking the symmetric-tracefree part of the 
relevant tensor. An alternative statement of the theorem written in the more natural language 
of space spinors will be given in the main text. It is a non-obvious result — to be proved in the 
main text — that the condition on the tensor Cij is actually a condition on the freely specifiable 
tensor $i,-. 



Overview of the article 

The article is structured as follows. In section [21 in order to fix notation and conventions some 
relevant preliminaries concerning the construction of spacetimes from an initial value problem are 
reviewed. These include the construction of asymptotically Euclidean solutions to the vacuum 
Einstein constraint equations from free data on a compact manifold. There is also a brief review 
of the solutions to the momentum constraint on flat space and also in the more general case where 
the conformal metric is not flat. Of particular relevance is subsection 12.41 where the assumptions 
made on the class of initial data being considered are spelled out. 

Section [3] contains a brief discussion of the construction of the bundle manifold Ca which 
provides a convenient alternative representation of the asymptotic region of the initial manifold: 
one in which the point at infinity i is blown up to a set which is topologically a 2-sphere. In 
particular subsections 13.31 and 13.51 discuss aspects of normal expansions at infinity. The ideas of 
these subsections are essential for the sequel. The techniques in section [3] are slight adaptations 
of the original constructions given in the seminal reference |11| — there is a more recent detailed 
discussion in [12j . The purpose of this section is to introduce necessary notation and to serve as 
a quick reference for the reader. 

Section m discusses the solutions of the momentum constraint of the bundle manifold Ca- This 
section builds from an analysis given in references [20] and [H]. However, the particular level 
of detail required for our analysis, in particular with regards to the non-flat case, is not to be 
available elsewhere in the literature. 

Section [5] builds on the results of the previous section to provide a detailed analysis of the 
normal expansions of the Weyl spinor near infinity. The main results of this section are presented 
in theorems [3] and m These theorems can be regarded as the generalisation of theorem 4.1 in [TT] . 
for initial data sets with a non-vanishing second fundamental form. 

Section [6] gives a brief overview of the ideas behind the so-called regular finite initial value 
problem at spatial infinity. Again, the discussion is kept to a minimum and has the purpose of 
introducing notation and ideas which will be of relevance in the sequel. In particular, subsection 
16.31 contains the crucial result of [11] stating that the solution to the conformal propagation 
equations will generically develop logarithmic singularities at the sets where null infinity touches 
spatial infinity. How to eliminate such singularities by means of conditions under the class of 
initial data under consideration is the concern of the rest of the article. 

Section [7| provides the main result of the article, theorem [G] which recasts the technical 
regularity conditions in terms of conditions on the Coton tensor of the conformal metric and on 
what is essentially the curl of the second fundamental form. This theorem brings together the 
discussions of section [S] and [Sj 

Section [8] discusses the possibility of generalisations of the main result to other classes of 
data. A particularly desirable generalisation would be to of a class of data including stationary 
solutions. 

Finally section [9] discusses the connections of the conditions obtained in the main result to 
similar — but crucially not identical — conditions which have been obtained from an analysis of 
purely radiative spacetimes — see [9] and also the concluding remarks in [TT] . 

The article contains two appendices. The first one contains a number of spinorial identities 
which are used through out the main text. The second appendix contains a discussion of the 
asymptotic expansions of the solutions of the vacuum Einstein constraint equations for the class 
of data under consideration. The appendix [B] briefly reviews some relevant results obtained in [5] 
and extends some of the analysis therein as the data required to render our analysis non-trivial 
turns out to be slightly more general than the one considered in the aforementioned reference. 

One of the main challenges of the analysis presented in this article is to bring together results 
and ideas obtained in different frameworks and cast them in a common language which allows, 
in turn, to obtain new results and hopefully also valuable new insights into the structure of the 
gravitational field near spatial infinity. 



2 Preliminaries 

This article is concerned with properties of asymptotically flat spacetimes {Ai , g) solving the 
Einstein vacuum field equations 

^M- = 0. (1) 

The metric g^j/ will be assumed to have signature (+, — , — , — ) and fj,,u, . . . are spacetime indices 
taking the values 0, . . . , 3 . The spacetime {A4,g) will be thought of as the development of some 
initial data prescribed on an asymptotically Euclidean Cauchy hypersurface S. The data on S are 
given in terms of a 3- metric hij of signature (—,—.—) and a symmetric tensor field Xij representing 
the second fundamental form induced by ^^^ on S. As mentioned in the introductions i,j, fc, . . . 
will be spatial tensorial indices taking the values 1,2,3. The Einstein vacuum field equations 
imply the constraint equations 

^'^r-x\+r'x^,=0, (2) 

D'X^, - D,r. = 0, (3) 

where D denotes the Levi-Civita connection and *^'^^f the Ricci scalar of the metric hij. The initial 
data {S, hij, Xij) will be assumed to be asymptotically Euclidean — for simplicity the situation 
with only one asymptotically flat end will be considered. In the asymptotically flat end it will be 
assumed that coordinates {y*} can be introduced such that 

with |yp = (y^)^ + (y^)^ + (j/^)^ and TO a constant — the ADM mass of 5. In addition to these 
asymptotic flatness requirements, it will be assumed that there is a 3-dimensional, orientable, 
smooth compact manifold (5, ft.), a point i G 5, a diffeomorphism $ : 5\ {i} — > S and a function 
n £ C^{S) n C°°(5) with the properties 

n{i)^0, Djn{i)^Q, DjDkn{i) = -2hjk{i), (4a) 

17 > on 5 \ {i}, (4b) 

hij = Vi^^Sij- (4c) 

The last condition shall be, sloppily, written as hij = fl^hij — that is, S \ {i} will be identifled 
with S. Under these assumptions (5, h) will be said to be asymptotically Euclidean and regular. 
Suitable punctured neighbourhoods of the point i will be mapped into the asymptotic end of S. 
It should be clear from the context whether i denotes a point or a tensorial index. 

2.1 The constraint equations 

In order to discuss the asymptotic properties of the solutions of the vacuum Einstein equations 
([T|) , the latter will be rewritten in terms of a suitable conformal factor f2 and a conformal metric 
g/iu such that 

gfiv = ' ' ff/^i/- 

The first and second fundamental forms determined by the metrics g^^i/ and g^^i/ on S are related 
to each other via 

h^j = fl'^hi-j, Xi] = ^{Xij + '^hij), 

where S denotes the derivative of fl in the direction of the future directed g-unit normal of S. If 
X = h'^^Xiji X = h^''Xijj one has that 

nx^x + 3S. 

In terms of the fields Q, hij and Xijj the constraint equations take the form 

2VlD,D'Vl - 3D^nD'n + -n^^^\ - SY."^ - -n^ (x^ - XvX'^) + 2^Sx = O, (5) 

n^D\fl-^Xtj) - ^ {DjX - ^Vt^^DjT) = 0, (6) 



where D denotes the Levi-Civita connection and ('^V the Ricci scalar of the metric h. In the 
sequel it shall be assumed that 

S = 0, x = on 5. (7) 

That is, the hypersurface S will be assumed to be maximal with respect to both /i.y and hij. 

2.1.1 The conformal constraint equations 

We now briefly review the conformal constraints implied by the metric conformal field equations 
on a hypersurface S — see e.g. [7] . The latter are written in terms of the fields g^^^ , Q and 

as 

2f7S' - V^f^Vf^ = 0, 

Va V - ^pLx. = V^nW^'^^p, 

where V^ denotes the Levi-Civita covariant derivative of /i^^, and R^j^vXp, R^iv, R are the as- 
sociated Riemann tensor, Ricci tensor and Ricci scalar of gfj^^, while C^uXp denotes its Weyl 
tensor. The tensor Wf^^Xp will be called the rescaled Weyl tensor. It will be of crucial relevance 
in the present article. The constraints implied by these equations on S will be discussed using 
a g-orthonormal frame field {ea}, a = 0,1,2,3 such that bq corresponds to the g-normal to S. 
Furthermore, let V^ = '^a^ p.- In what follows, the (frame) indices a,b,c,.. . will be assumed to 
take the values 1, 2, 3. Set J^ — Vq^I and W*,^^^ — ^Wpiyape"'^ . One further defines 

Wabcd = Wabcd: Wab = WaObO, W*^ = W^obO' '^abc = WaObc, 



satisfying 



= -2 {ha[cWd]b + h^dWc]a) , "'adS'^foc = ^abc, W^^ = -^Wabcf^J""^ 



Wabcd 

Wab = Wta, W°'^ = 0, Wlt,=wl^, U;*° == 0, 

Wabc = -Wacb, w"- ^^ = 0, W[abc] = 0, 

where the indices are moved with hab ~ ~Sab, and tabc is the totally antisymmetric 3-dimensional 
tensor with 6123 — 1. The tensors Wab and w*^ represent, respectively, the n-electric and n- 

magnetic parts of W^ l^,; on S. Using the Gauss' and Codazzi equations on S one obtains the 
following interior equations (the conformal constraint equations): 

2nS-J:^ - DanDTl = 0, (8a 

DaDbfl = -SXafc - ^Lab + Shab-, (8b 

Dj: = Xa'Dc^ - ilLa, (8c 

DaS - -D'nLta - Sia, (8d 

DaLbc - DbLac = D^flWecab - ^Wcab - XacLb + XbcLa, (8e 

DaLb - DbLa = D^flWeab + Xa ^bc ^ Xb^^ac, (8f 

D^Wcab = X^aWbc - X^b^ac, (8g 

D^Wab = X^'Wabc (8h 



These equations can be read as conformal constraints for the fields Q, S, S, hah, Xat, La, Lab, 
Wab, Wab- Alternatively, — and this is the point of view which shall be adopted here — if one has 
physical data, {hij,Xij) solving the vacuum constraint equations and a conformal factor i7, and 
moreover, a choice of E and R has been made, then one can use the equations to calculate the 
initial data for the conformal field equations. 

2.2 Solving the constraint equations 

Under the assumptions ((7|) the conformal constraint equations take the form 

U _ 1(3)^^ ^ = \x^x''^, with d = ^-^'\ (9) 

D' {n-\,,) = 0, (10) 

with A = Dj^D^ . In what follows it shall be assumed that the above equations are solved using 
an adaptation of the so-called conformal method. Namely: 

1. choose a smooth, negative definite metric hij on a 3-dimensional, oricntable, smooth com- 
pact manifold S and pick a point i (z S. Set 5 = 5 \ {i}. 

2. Find a smooth, symmetric tensor field V'lj on S which is trace-free with respect to hij and 
satisfies 

D"^^J = 0. (11) 

The tensor ipij can be obtained by means of a York-splitting. Given a smooth, symmetric. 



trace-free tensor $y on S, set 



2 J. 

i/jij = DiVj + D-jVi - -hijDkV + $,j. 



iCv),j + $, 



u 



where Vi is some 1-form on S and {Cv)ij is the conformal Killing operator of hij. Given the 
above Ansatz, then equation ()lip implies the following elliptic equation for Vi'. 

Av, + ^DjDkv'' + rjkv'' = -D'^^^,, (12) 

which under suitable conditions can be solved. 
3. Setting Xij — i? '^V'ij in equation ^ one obtains the Licnerowicz equation 

(^A - i(3)^^ ^ = iv'.,^^^l?-^ (13) 

that is, an elliptic equation for i!). 

The fields hij, il = -i?"^ and Xij = ^'^ipij so constructed render a solution to the conformal 
constraints ([5]) and ^. It is important to recall that if is a positive scalar field on S then the 
transitions 

hij ^ (jy^h.j, tp.j —> (jy^'^ij^ij, Vi-^cjP'VL, Xij ^ <f>'^Xi3, (14) 

provide another solution to the conformal constraint with the same physical data hij, Xij- This 
freedom in the conformal gauge will be used in the sequel to obtain simplifications in the analysis 
of asymptotic expansions. 

Given a 3-metric hij consider a suitable a > such that Ba = Ba{i) is a strictly convex normal 
neighbourhood of i and let x* be normal coordinates with origin at i based on an /i-orthonormal 
frame Ca- The asymptotic condition (j4ap implies 



|a;|i? ^-1 as a; ^- 0. 



Accordingly, one has that 

Xy = C'(l), Vjj = f -TjI j as a; -^ 0, 
where the following transformation rules have been taken into account: 

Under suitable conditions, the conformal factor -d admits, in a neighbourhood Ba with a > 0, 
the decomposition 

^=^ + W, (15) 

where 

with S(i) the Dirac's delta distribution with support on i and 

The function C/ describes the local geometry in Sa and can be calculated by means of Hadamard's 
parametrix construction. If hij = —^ij — the conformally flat case — then U = 1. The function 
W contains global information about {S,hij,Xij)- 

2.3 The momentum constraint in flat space 

The analysis of the solutions of the momentum constraint in flat space given in [5] is now briefly 
recalled. In section |4] it will be recast in a form — in terms of space spinors — which is convenient 
for the application discussed in this article. 

In this section assume that S is flat in at least a neighbourhood Ba of i. Let x'^ be a Cartesian 
coordinate system with origin at z. In these coordinates one has that the metric is given by 
hij — —Sij. We shall often write r — \x\. Let n^ = x^/\x\ and complement it with complex 
vectors m% m* to form a basis of the tangent bundle of R^, by requiring that 

rriim^ — rnifrf — riirn^ — nirrf ~ 0, mifn^ = — 1. 

Recall that in this construction one has freedom of performing rotations m^ i— > errii with 9 
independent of r. The metric hij can be written in the form 

hij = riiUj + niimj + mimj, 

while an arbitrary trace-free tensor ■0^ can be written as 

ipij — — (^{SniUj — Sij) + V2r]n(^ifnj^ + v^nf^^rrij) + JLmimj + firnirrlj j , (16) 

with 



? 



—ri/jijU^n-', rj — V2r ipijU^m-' , fx = r ipijm^m- 



^mJ 



2 

In the previous expression ^ is real, while rj and /x are complex functions of spin weight 1 and 2 
respectively. Let A be an arbitrary complex C°° function on Ba- Let A = 3^ A, where S denotes 
the ei/i-operator —see e.g. [UlIlT]. Let A^^) = Re(A) and A^-^) = ilm(A). Let 

^'j ^ A [-P'n^ - pin' - {S'^ - 5n'n^)P''nk) , (17a) 

rj = A {n'e'''^J,ni + n^e'^^' J^n,) , (17b) 

^J = 4 (3nVi^' - S'^) , (17c) 

i;^ EE ^ {QV + Q^n' - {S^^ - n'n^)Q''nk) , (17d) 



where P' (the linear momentum), J' (the angular momentum), Q^ (conformal momentum) and 
A are real constant vectors, respectively, scalars. Furthermore, let 

^^tx^^'\ (18a) 

r] ^ ~2rdrW^^ +W^\ (18b) 

H = 2rdr{rdrX^^'>) - 2X^^^ + dEX^^'> - rdrX^'\ (18c) 

and denote by V'ij[A] the tensor obtained by substituting (|18ap - (|18cp into (fT6)l . Then 

tP,j = ^ij [P] + l/)y [J] + Viy [A] + i^ij [Q] + tP,J [A] , (19) 

satisfies the flat space momentum constraint 

d% = 0. (20) 

Conversely, any smooth solution to equation (|20p is of the form P^ — cfr. theorem 14 in [5]. 

In the present article the complex function A which is smooth on Ba \ {i} will be taken to be 
of the form 

A = A(i' + iA(2)^ (21) 

r 

with A(i), A(2) e C°°(Sq). In analogy to equation ^ let also 

r 

Note that the factor 1/r makes the function A non-smooth at i. The functions A'^-* and A*-^' are 
calculated using the 3 operator. Accordingly, they are non-smooth. 

2.4 Assumptions on the freely specifiable data 

In order to obtain a non-trivial outcome from the analysis to be discussed in the sequel, one 
has to choose a suitable class of initial data. If the initial data is to be constructed following 
the procedure discussed in section 12.21 then the freely specifiable data is given in terms of a 
negative-definite 3- metric hij and an ft,-trace free symmetric tensor $ij. 

2.5 The conformal metric 

For ease of the presentation, the results presented here will be restricted to the case of conformal 
metrics hij which are analytic in a neighbourhood, Ba, of spatial infinity. This assumption is 
non-essential and with long recursive arguments it should be possible to extend all the results 
presented here to the smooth — i.e. C°° — setting. The analyticity of hij is explicitly used in the 
proof of part (ii) of theorem [3] and of part (i) of theorem |4l where the expansions of a certain 
tensor associated to the massless part of the Weyl tensor are calculated. 

It should be mentioned that in some senses the assumption of analyticity of the conformal 
metric on Ba is as restrictive as that of smoothness. The reason behind is that generic stationary 
spacetimes do not have a smooth conformal metric. Instead, as shown in [3], the metric is of the 
form 



-ij — 'Hj < ' '"ij 



.(1) .„A ^(2) 



with hi-' and h^-' analytic tensors. Due to the presence of the term r"^, the metric in this case 
will be C'^'°'{Ba), < a < 1. The regularity of this conformal metric is conjectured to be optimal 
for strictly stationary data — i.e. stationary data which is not static. 

Throughout, it will be assumed that in normal coordinates centred at i one has that 

h,^=-S,j+h,j, h"^ =-5"^ +h}\ 



with hij e C'^{Ba), /i*-' G C'^{Ba)- It is recalled that the normal coordinates satisfy 

x^h^j = 0, Xj^*-' = 0. 

The latter relations imply that 

Xix'V\j = 0. 

The analysis of the solutions to the momentum constraint performed in the appendix |B] requires 
the technical assumption 

S'^h.j ^S,jh'^ ^0. (22) 

The conformal gauge freedom p^ can be used to obtain a representative in the conformal 
class of hij which is in the so-called en (conformal normal) -gauge — see jll| . Working with a 
metric in the cn-gauge renders a number of useful simplifications in the present discussion. In 
order to define the cn-gauge, consider on Ba solutions (x^ (t) , bi{t)) of the eonformal geodesie 
equations 

x'^Dkx' = -2{bkx'')i' + {hkix''x')h"''b„,, (23a) 

x'^Dkb, = (6fci'^)6, - ^{h'''bkbi)h^,x"' + ^,fei^ (23b) 



with initial conditions 



x{0) = i, hijx'x^ = -1, 6,(0) = 0, 



where x^{t) denotes a curve in Ba through i and bi{t) is a 1-form along that curve. If neigh- 
bourhood Ba of i is taken to be small enough, then there is a unique conformal rescaling of the 
form given by (|14p such that the conformal metric is analytic, keeps the metric and connection 
unchanged at i, and such that if the solution to the conformal geodesic equations (|23ap and (|23bp 
is written in terms of the rescaled metric one has 

b,x' ^ 0, on Ba. (24) 

A metric in the conformal gauge of hij satisfying (j24p will be said to be in the cn-gauge. From 
here on, it will always be assumed that the metric hij is in the cn-gauge. 

An important consequence of working in the cn-gauge is that both the Ricci scalar and Ricci 
tensor of hij vanish at i. Accordingly, 

kj = 0{r^), duK, - 0{r''). 

In addition, it will be proved — see lemma [5] — that if hij is in the cn-gauge then hij — r^hij so 
that 

hij ^ ~b,j + r^K.j , K^ = -(5*^ + r^/i^J , 

where hij — 0(r) is analytic. Note that because of the condition ((22)) one has that 8'^^hij = 0. 

2.5.1 Free data for the second fundamental form 

The /i-trace free symmetric tensor $ij , containing the free data for the second fundamental form 
will be constructed out of the tensor ipij , solution of the flat space momentum constraint given 
in equation (fT9)) . Thus, let 

$,j =ipij - -hijh'^^ipki. 

Thus, $ij is specified by prescribing the constant vectors P% J\ Q*, the constant scalar A, and 
the function A. In order to obtain a vector w* and a scalar 'd solving the equation (J13p admitting 
an expansion in powers of r near i one has to set P^ — and consider a A of the form given in 
equation (|2ip . A more extended discussion of these issues is given in [5]. However, it turns out 
that if one admits a contribution to the second fundamental form of the type ipij with A given by 
equation pTjl . the solutions to the elliptic equation (fT2|) will not have an asymptotic expansion 
in Ba consisting purely of powers of r. In order for this to be the case, Im(A'^^^) and Re(A^^-') 



have to be related to the conformal metric hij in a particular way. This issue is discussed in more 
detail in the appendix. Throughout the main body of the article, it will be assumed that such 
conditions are satisfied, and accordingly, the solutions w* of equation p^ have expansions in Ba 
purely in powers of r — see below for the precise details. Examples of classes of data where these 
are valid are axially symmetric data and conformally flat data — in which case f ' — 0. 

2.5.2 Some consequences of the assumptions on the freely specifiable data 

In order to make precise the idea that a given function or tensor field over S admits an expansion 
around i in terms of powers of r, introduce the following function spaces: 

E°°{Ba) = {/ = /i + r/2 I /i, /2 e C°<-\Ba)} , 

QooiBa) = {v' e C°^{Ba,R') I xy = r^v, V e C^{Ba)} . 

Theorem 15 in [5] states that if the function A determining the higher multipoles of tpij is 
of the form given by equation ^ —i.e. Xr e E°°{Ba) and P^ = 0, then r^-ijijip-^i e E°°{Ba), 
from where it follows from their main result — theorem 1 — that in the conformally flat case 
{hij = —Sij) one has that 

^^1/r + W, W€E°°iBa). 

In the non-conformally flat case, if hij and A are such that the solutions, w*, of equation (fT2|) are 
of the form 

for some integer s with v\ e QooiBa), v^ e C°°(;Ba), and such that r^^yV'*'' ^ E°°{Ba), with 

7/.y = 4„-[A J,Q, A(i\ At^Vr] + {Chv),,[A,J,Q] + (£;,t;)„[AW] + (£,,«)„ [A^^)], 
then theorem 1 in [5] renders 

r 
with U S C^{Ba) and W G E°°{Ba)- Furthermore, because of the use of the cn-gauge, U = 
1 + 0{r'^). The discussion in appendix [B] renders a detailed description of the structure of the 

vectorsw*[A, J,Q,A(i),A(2)/r]. 

3 The Manifold Ca 

In [11] a representation of the region of spacetime close to null infinity and spatial infinity has 
been introduced. The standard representation of this region of spacetime depicts i" as a point. In 
contrast, the representation introduced in |11J depicts spatial infinity as a cylinder — the cylinder 
at spatial infinity. The technical and practical grounds for introducing this description have been 
discussed at length in that seminal reference. The original construction in [llj was carried out 
for the class of time symmetric metrics with analytic conformal metric hij. However, as discussed 
in [201 m] , the construction can be adapted to settings without a vanishing second fundamental 
form. The purpose of this section is, primarily, to introduce notation that will be used in the 
sequel and to provide enough background material to follow the discussion in the sequel. In any 
case, the reader is referred to [Til IE] for a thorough discussion of the details. 



3.1 The construction of the manifold 

Starting on the initial hypersurface S, Friedrich's construction makes use of a blow up of the 
point i G 5 to the 2-sphere § . This blow up requires the introduction of a particular bundle 
of spin- frames over Ba- In what follows a space spinor formalism analogous to a tensorial 3+1 
decomposition will be used to this end. Consider the (unphysical, conformally rescaled) spacetime 
{A4,g^,y) obtained as the development of the initial data set {S, hij, Xij)- Let SL{S) be the set of 
spin dyads 5 = {Sa}a=o,i = {oa, m} on S which are normalised with respect to the alternating 
spinor eab in such a way that eoi = 1- 
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The set SL{S) has a natural bundle structure where S is the base space, and its structure 
group is given by 

5i(2,C) = {t^s e GL{2,C) I eAct^st'^D = ^bd}, 
acting on SL{S) hy S i-^ 6 ■ t ^ {'5a^"^s}b=o,i- Now, let r = V^eo, where cq is the future 5-unit 
normal of S and 

TAA' = 9(t, Sa'Sa') = e/e^,° + e/e^^} 
is its spinorial counterpart — that is, r — r'^ea — '^aa'''' ^a where c'^a' denote the Infeld-van 
der Waerden symbols and {ca}, a = 0, ... ,3 is an orthonormal frame. The spinor taa' enables 
the introduction of space-spinors — sometimes also called SU{2) spinors, see [TJinillS]. It defines 
a sub-bundle SU{S) of SL{S) with structure group 

SU{2,C) = {t^B e SL{2,C) I TAA't^Bt^'s' = rsB'}, 

and projection n onto S. The spinor t"^^ allows to introduce spatial van der Waerden symbols 
via 

AB {A B)A' a A' a ■ i o Q 

^a ^^a A'^ ' ^AB=^(S ^A)A': i^ 1,2,6. 

The latter satisfy 

7 AB cab 7 

riab = O-aABO-f, , —OabCTAB'^CD = —<^A{C^D)B — IT-ABCD, 

with /i of, = h{ea,eb) — —Sab- The bundle 5^/(5) can be endowed with a 5u(2, C)-valued connection 
form u! B compatible with the metric hij and 1-form <t , the solder form of SU{S). The solder 
form satisfies by construction 

h = h^jdx' (g) dx^ = hABCDO^^ ® cr^", (26) 

where cr^^ = af^dx"^ — note that the crf^ arc not the spatial Infcld-van der Waerden symbols. 

Now, given a spinorial dyad S G SU{S) one can define an associated vector frame via Ca — 
ea(^) = '^t^^ATg^ 5b', 0, = 1,2,3. We shall restrict our attention to dyads related to frames 
{sj}j=o,-- ,3 on Ba such that 63 is tangent to the /i-geodesics starting at i. Let H denote the 
horizontal vector field on SU{S) projecting to the radial vector 63. 

The fibre 7r^^(i) C SU{S) (the fibre "over" i) can be parametrised by choosing a fixed dyad 
5* and then letting the group SU{2,C) act on it. Let {—a, a) 3 p t-^ 5{p,t^g) G SU{S) be the 
integral curve to the vector H satisfying (5(0, t"^^) = 5{t^ g) G 7r^^(i). With this notation one 
defines the set 

Ca = {S{p,t^s) e SU{Ba) I \p\ < a, t^B G SU{2,C)}, 
which is a smooth submanifold of SU{S) diffeomorphic to {—a, a) x SU{2,C). The vector field 
H is such that its integral curves through the fibre Tr~'^{i) project onto the geodesies through i. 
From here it follows that the projection map tt of the bundle SU{S) maps Ca into Ba- 

Let, in the sequel, X° = 7r~^(i) = {p = 0} denote the fibre over i. It can be seen that 
X" « SU{2, C). On the other hand, for p E Ba\ {i} it turns out that tt~^{p) consists of an orbit 
of C/(l) for which p = \x{p)\, and another for which p = —\x{p)\, where x*(p) denote normal 
coordinates of the point p. In order to understand better the structure of the manifold Ca it is 
useful to quotient out the effect of f/(l). It turns out that I''/f/(l) ~ S^. Hence, one has an 
extension of the physical manifold S by blowing up the point « to § . 

The manifold Ca inherits a number of structures from SU{S). In particular, the solder and 
connection forms can be pulled back to smooth 1-forms on Ca- These shall be again denoted by 
a and uj g ■ They satisfy the structure equations relating them to the so-called curvature form 
determined by the curvature spinor 

rABCDEF = I 7;^ABCE — J^rhABCE HDF + I -^^ABDF — JT^rtlABDF j^CE, 

where sabcd = S(^abcd) is the spinorial counterpart of the tracefree part of the Ricci tensor of 
hij and ^^'r its Ricci scalar. These satisfy the 3-dimensional Bianchi identity 



D^^ SABCD = ^DcD^'^^r. 
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3.2 Calculus on Ca 

In the sequel t^ and p will be used as coordinates on Ca- Consequently, one has that H — dp. 
Vector fields relative to the SU{2, C)-dependent part of the coordinates are obtained by looking 
at the basis of the (3-dimensional) Lie algebra su(2, C) given by 

Oi\ 1/0 -1\ 1/iO 



"1 = 2 1, i ; ' "2 ^ 2 1, 1 ; ' "3 = 2 V -i 

In particular, the vector u^ is the generator of U{1). Denote by Zi, i = 1, 2, 3 the Killing vectors 
generated on SU{S) by Ui and the action of SU{2, C). The vectors Zi are tangent to X'-^ On I^ 
one sets 

X+ = -{Z2 + iZi), X_ = -(Z2-iZi), X = -2iZ3, (27) 

and extend these vector fields to the rest of Ca by demanding them to commute with H — dp. 
For latter use it is noted that 

[X,X+] - 2X+, [X,X^] = -2X_, [X+,X_] = -X. 

The vector fields are complex conjugates of each other in the sense that for a given real-valued 
function W, X-W = X+W. More importantly, it can be seen that for p E Ba\{i} the projections 
of the fields H, X± span the tangent space at p. 

A frame cab — C(ab) dual to the solder forms cr*^^ is defined so that it does not pick 
components along the fibres — i.e. along the direction of X . These requirements imply 



cr 



, ccd) = h^^cD^ ccD = ccijdp + c^jjX+ + c^^X^, (28) 

where (•, •) denotes the action of a 1-from on a vector. Let a^ and a be 1-forms on Ca annihilating 
the vector fields dr, dp and having with X± the non- vanishing pairings 

{a+,X+) = (a",X_) = {a,X) = 1. 

Furthermore, from the properties of the solder form cr one finds that 

Cab = XAB, c\g = -ZAB + c\b, c^b = -Vab + c^b^ (29) 



with constant spinors xab, VAB and zab given by 

1 

72^"^"' •'"" V2 



XAB = \'20(^A>'B), VAB ^ -—7=I^Al^B, ZAB ^ —7=OAOb, 



and analytic spinor fields satisfying 

cSb=0(p), c^i=0, a=l,±. 
Accordingly, one can write 

Cab = CyyAB + C^ZAB- 

Furthermore, using the structure equations it can be shown that in fact c\q = 0. By virtue of 
the relations (|28p and (I^S)) . the solder forms a^^ descend to forms riidx*, m^dx*, mMx' spanning 
the tangent space of the points of Ba with non- vanishing pairings 

•n}ni = — 1, mim^ = — 1, 

and such that n* = x^ /r. Note that n*, m^ and nf are not smooth functions with respect to the 
coordinates a;*. 

The connection coefficients are defined by contracting the connection form Co b with the frame 
CAB- In general, one writes 

ICD B = {<^ B^Ccd) = -ICD B +lCD B: 
lABCD = -^i'^ACXBD + ^BDXAc)- 
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The spinor 7^bc£> denotes the singular part of the connection coefhcients. The regular part 
of the connection can be related to the frame coefficients cab via commutator equations. The 
smooth part ^abcd of the connection coefficients can be seen to satisfy 

711CD = 0, JABCD = 0{p). 

Furthermore, from an analysis of the structure equations one obtains that 71100 = ~7ooii- Let / 
be a smooth function on Ca 

Dab! ^cabU)- 

Similarly, let ^lab represent both a smooth spinor field on Ca- Then the covariant derivative of 
^1AB is given by 

DABfJ-CD = CABifJ-Co) - IaB C^^ED - 1 AB D^^CE■ 

Analogous formulae hold for higher valence spinors. 

3.3 Normal expansions at i 

In [11] a certain type of expansions of analytic fields near i has been discussed. Although the 
conformal metric will be assumed to be analytic on Ba, most of the other fields — in particular 
objects derived from the second fundamental form — will at most be smooth. Thus, the ideas of 
[TT] have to be adapted to the smooth setting. This can be readily done. 
If / e C^iBa) then one has that 



/^E/*- 



■i^X 



fe>0 



which is to be interpreted as 



/ = E/,,...,,x'^---x''=+/,,, 



fe=0 

with Jr e C^{Ba), Ib. = oir"'), for all m > 0. The term EI'Lo /n «. 2:'^ • • • x'% with /,,...,, 
constant vectors, is the Taylor polynomial of degree m of /. 

Now, suppose that yi''"- is a smooth tensorial field of rank(r, s) on Ba with components 

rp*ai---ar^ ^ with respect to the frame e* on which the normal coordinates a;* are based. Let 
V = x->dj be the radial vector which is tangent to geodesies through i and satisfying ViV"^ = —1. 
Let also ri^ be defined by V"^ = |a;|n'. By construction one has that V'^Dke'a = 0, thus following 
the procedure described in section 3.3 of [TT] one obtains 

p=0 ^' 

with q e Ba, and iT*'''-\...jTz G C^iBa) and iT*'''-\^,„Jn = o^r^). The first term in 
the right-hand side of expression pO|) will the called the analytic part of x*°'^"'°"'^^ ...j^ . 

An analogous expansion can be obtained for smooth analytic spinor fields. Suppose ^AiBi---AiBi 
denotes the components of a smooth even rank spinorial field with respect to the spin frame 5\ 
associated to e*. In analogy to expression (|30p one can introduce the expansion 

771 . 

Cmb.-MbM) - E -M\'n^'''^ ■ ■■n^'^^^Dc^n, ■ ■ ■ Dc,d,Cmb,...a,b, + {Cmb.-mb,)^ ' (31) 
p=0 ^' 

with n^^ ~ v^^{q), q £ Ba and the derivatives of the spinor field are evaluated at the point i, 
and {Ca,Bi-AibM e C°°{Ba), {C1,b,-AiBi)'R- = o(r™)- Again, the first term of the right-hand 
side of equation (j3ip will be referred to as the analytic part. 

An analysis of the decomposition in terms of irreducible spinors of the summands in (|3ip will 
be important in the sequel. The derivatives in the expression (|3ip can be replaced by symmetrised 
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derivatives (in the pair of indices CkDk) ^^ they are contracted with the same spinor n . In 
these symmetrised derivatives, the contraction of indices Cj j Ck with j ^ k renders an expression 
antisymmetric in the indices d and d^. ■ Accordingly, the decomposition of every sumand in (j3ip 
with p > 2 into symmetric, irreducible parts — with respect to the indices c^Dp-CiDi — renders 
an expansion 

„C,D, . ..^C,D,jj^^^^ . ..Dc,D,Ca,B,-A,B, 

' {ip,0;CpDp-CiDiAiBi-AiBi + ^p,l:CpDp---C3D3AiBi---A,Bi^CiDiC2D2 + 

where the dots indicate terms at least quadratic in hABCD- In particular one has that 

^p,0;CpDp---CiDiAiBi---A,B, — D(CpDp ' ' ' -DciI3i)CaiSi ■■■A,Bi i 

Cp,0-CpDp-CiDtAtBt-AiBi " 2n - 1 ^ D^^CpDp ' ' ' D\ef\ ' ' ' D ' ' ' Dc^D3)CAtBt-AiBii 

^ l<k<h<p 

with Def and D^^ assumed to be, respectively, in the kth and hth position. More generally, 
one has the symmetry 

^pA;CpDp---C2i+iD2, + iAiBi---AiBi = ^p,i;{CpDp---C2.+iD2i+i){AiBi---AiBi)^ 

with i = 0, . . .p/2 if p is even and i — 0, . . . , {p — l)/2 if p is odd. The complete decomposition 
of the coefficients in ([31]) renders 

^* AiBi---A,B, 

^p,i;CpDp---C2i+iD2i+i 

_ tO AiBf-AiBi f.1 (AiBf-A, Bi) 

^p,i;CpDp--C2,+iD2i+i ~^ ^p,i;(CpDp---C2i+i D2i+i) 

,f2 (AiBi--B,_i A, Bi) ,„„^ 



with 



ip,i;CpDp---C2i+iD2,+iAiBi---AiB, - ip,i;{CpDp---C2^+iD2. + iAiBi---A,B,)^ (33) 

^p,i;CpDp---C2i+iD2i+iAiBi---AiBi ~ ^p,i,lip,i-{CpDp---C2i+i\E\AiBi---Ai) ' (34) 



for some real coefficients cl ^ i- 



3.4 An orthonormal basis for functions on SU{2, 



The lift of the expansion (|3ip from Ba to Ca introduces in a natural way a class of functions 
associated with unitary representations of SU{2, C). Namely, given t^^ G SU{2, C), define 

/ \ 1/2 / X 1/2 

^0 o(* s) = Ij 

with j,k — 0, . . . ,m and to = 1, 2, 3, . . .. The expression (Ai • • • Am)k means that the indices are 
symmetrised and then k of them are set equal to 1, while the remaining ones are set to 0. Details 
about the properties of these functions can be found in [51 [TT]. The functions y/m + IT^\ form 
a complete orthonormal set in the Hilbert space L^{fi, SU{2, C)), where fi denotes the normahsed 
Haar measure on SU{2,C). In particular, any analytic complex- valued function / on SU(2,C) 
admits an expansion 

oo m m 
m=0 j=0 k=Q 

with complex coefficients fm,k,j- Under complex conjugation the functions transform as 



rn k V / m m — k' 
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The action of the differential operators (|27|) on the functions T^ 'l- is given by 



A function / is said to have spin weight s if 

Xf = 2s f. 
Such a function has a simplified expansion of the form 

oo m 

J ~ / J / , Jm.k-Lm m/2-s- 
Tn>\2s\ fc=0 

3.5 Normal expansions at X° 

In the sequel it will be necessary to be able to relate fields in Ca with fields in Ba- Crucially, one 
will require to be able to lift smooth fields defined on Ba to Ca- As in section p.3|) consider normal 
coordinates x* on Ba centred on i based on the orthonormal frame c* — cr^^c^g = a^^S'XS'^. In 
terms of p and t g on Ca and the normal coordinates x% the projection tt' has the local expression 

tt' : ip,t^g) ^ x\p,t^g) = V2pahDt''ot°i- 

This expression can be used to pull-back fields (lift) to Ca- In particular the pull-back of |a;| is p. 
The procedure to lift expansions of the type given by ((3T|) has been discussed in [TT] . Starting 
from the analytic part of ((3T|) one has to perform the following operations: 



(i) the transition Gi-A, ^ Cb^-b,^'^^, ' ' '^'^'a,; 

(ii) the replacement |x| -^ p and n^^ -^ \/2t ^t '^; 

(iii) decompose all the spinor- valued coefficients into sums of products of symmetric coefficients 
with e^s's. Contractions of e^s's with a pair of t'~^ ^ yields factors of or 1, while the 



D 

remaining expressions assume the form of expansions in terms of the functions T^ 



T7l / ■ 



Applying the previous procedure to the analytic part of the expansion (|3ip of a spinorial field 
S,AiBi---AiBi on Ba One obtains the expansion of the spinor-valued function ^AiBi---AiBi on Ca- 
Denote by ^j — £,{AiBi-AiBi) ■: < j < Z its essential components. The function S^j has spin 
weight s = I ^ j and a unique expansion of the form 

m 

(,^yj,.pP'' + i(,)n, (35) 



p=0 



for all m, with 



P+l 2q 

q—mayi{\l~j\,l—p} k—0 

and complex coefficients £,j^p-^2q,k- In particular one has that 
^ (2p + 2lV'^(2p + 2l\-"\ 

4i,p;29+2I,fc - (V2) 1^ ^ j Vp + jj ^iC^D,---Oc,D^kA,B^-ApB^)W- 

Hence, one has the symmetry 

Co,p;2p+2i,fc = £,2l,p;2p+2Lk 

which will play an important role in the sequel. Another lengthy, but straightforward calculation 
shows that 

^j,p;2p+2M = Kp,j^kD(^ApBp ' ' ' ^Ai\E\^ABC).i (*)' (2'^) 
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with Kpj^k a constant depending on p, j, k. If / is even, then S,Ai---Ai is associated to a real spatial 
tensor if and only if the following reality conditions hold: 

Cj = (~l)jf2;~j' ?J,p;2g,fc = (—1) ^2l~j,p;2q,2q-k- 

The discussion in the following sections will require to consider smooth spinorial fields £,AiBi---Ar.Br- 
with essential components ^j = K{AxBx---ArBr) i !i J ^ 2r of spin weight s — r — j with expan- 
sions which are more general than those given in equations (|35p and (|36p . Accordingly, they do 
not descend to smooth spinors on Ba- In this case, instead of (j36p . one considers the more general 
expression 

9(p) 2 

ij,P ~ 2^ / ,Kj,p;2q,k-L2a g^r+i' 
q=\r-j\ fc=0 

where in principle < |r — j| < q{p) < oo. In this case one will speak of an expansion of type q{p). 
If /, g have expansion types qip), q'{p), respectively, then the product fq will have expansion 
type maxo<j<p{q'(j') + q{p — j)}, while the sum f + g will have expansion type max{g(p), q'{p)}- 
In particular, pf will have expansion type q{p) — 1. 

3.6 Consequences of the cn-gauge for normal expansions 

In [11] a number of consequences for quantities derived from an analytic conformal metric hij in 
the cn-gauge has been deduced — see lemma 4.7 in the aforementioned reference. 

Lemma 1. In the cn-gauge one has that 

type(r)=p, r^O{p^), 

type(sABCD) = P + 1, SABCD = 0{p), 
type(7ABCD) = P, lABCD = 0{p^), 
type(c^B) = P, c^B = O(P^), 
type(C/-l)=p-l, U = l + 0{p^). 

In addition, from the discussion in section 12.5.21 one has that under the assumptions being 
made it follows that the function W in equation ((15]) satisfies W £ E°°{Ba)- The procedure of 
subsection 13.51 gives 



Lemma 2. 

Ttlj 

type(W^) = p, W=- + 0{p) 

Remark. It is worth noticing that by means of an adequate choice of the centre of mass it 
is possible to set W — m/2 -\- 0{p^) — see e.g. [20], although this fact will not be used here. 

An important consequence of lemma [H of crucial relevance in the analysis of solutions of the 
momentum constraint and in particular of equation (J12p is the following lemma already mentioned 
in section [ 



Lemma 3. In normal coordinates based around i an analytic metric in the cn-gauge is of the 
form 



h ■ — — f5- + r'^h- 



with hij — 0{r) 



Proof. As seen in section[21 the lift of the conformal metric hij to Ca can be expressed in terms 
of the soldering forms (j^^ as h ~ hABCD<y^^ ® a^^ . Now, using the explicit decomposition 



jAB ^ _2;ABj^ ^ (-2/9y^^ + <T^y^^ + aU^'^)a+ + (-2pz^^ + a^z^^ + a^y^^)c 
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one finds that 












-l- 




~ PCy Ct - C+) , 


-. 2p__ 




crl-~ 




2p\ .__+ 


' PCyCt -Cy) 



with 



/- 2 h 2 h -+ — 1 

? = p c^ S - p Cj^ cj - pc;^ - pc^ - 1. 



In the en-gauge c^g has expansion type p and hence a^^ has expansion type p—2. Furthermore, 
because in the cn-gauge c^^ — 0{p^) one concludes that cr^^ = 0(/5^). Accordingly one can 
write 

~AB 3~AB ~AB mt \ 

Now, it is noticed that 

h = hABCD {-x^'^Ap - 2p/^a+ - 2pz^''a-) ® (-x^^dp - 2py^''a+ - ^pz^'^'cT) , 

= -dp (g) d/9 - 2/9^(a+ (g) a~ + a" a^), 

on Ca descends to the flat metric (in Cartesian coordinates) 

h = —6ijdx^ (g) da;^, 

for the pull-back of the standard metric on the 2-sphere is given by 2(a+ ^a^ + oT ® a+). Thus 
the 1-forms —x^^dp^ —2py^^a^ and ~2pz^^ a'^ on Ca descend, respectively, to 1-forms on So, 
riida;', rhidx^ and rfiidx'^ with n' = x'^/r and such that 

h = (niUj + rhifhj + ifijTTiij dx* (g da;-'. 

Similarly, the 1-forms a^^a^ and (T_a^ on Cq descend to the 1-forms r^mjda;-' and r^rhjdx^ 
on Ba, with rhj and m^ analytic. Thus, one has that 

h = (nifij + rkirkj + rkjrki) dx^ (E) dx-' + r [rhiifij + niirrij + rhjihi + rhjfhij dx^ (E) dx-' 
+r* {ihifiij + rhjrhi) dx* ® dx-'. 

Now, rhi — 0(1) is not a smooth function of the normal coordinates, however, by assumption hij 
is smooth, and consequently the combination 

ihirhj + rhifhj + fhjfhi + rhjifii, 

must be analytic on Ba- Hence, the metric is of the required form. D 

4 The solutions to the momentum constraint on Ca 

In this section the lifts of the solutions to the momentum constraint introduced in section [2?3l will 
be analysed. The procedure described here has been previously implemented in [20ll22j. 

4.1 The conformally flat case 



In section 12.31 the solutions of the flat space momentum constraint have been discussed by intro- 
ducing a certain frame {ri\ m^, m*}. This frame is related to an orthonormal frame by means of 
the relations 

Inverting these relations one can readily rewrite {n*, to*, m*} in terms of the orthonormal frame 
{eo}a=i,2,3. This leads to a direct transcription into spinorial objects. Let {ua, ma, rha] denote 
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the components of the frame {n*, to% m*} with respect to {e^}, a — 1,2,3. Using the spatial 
Infeld symbols one obtains to the following transcription rules: 

na -^ XAB, rua ^ V2zab, ma -> V2yAB- 
Thus, formula dH]) lifts to 

" "IJABCD = <;\o-i.A. 

+2'Jl2iyABycD) + 2h2{zabzcd)- (38) 



P'^i'ABCD = ^{iXABXCD + hABCo) + V^rj^ixABVCD + XcdVAb) + V2l]iixABZCD + XcdZAb) 



Alternatively, one could rewrite the previous formula in terms of the totally symmetric spinors 
^ABCD^ * = 0, 1, . . . , 4 as: 

P^i'ABCD = IpO^ABCD + '^I'^ABCD + '^2^^300 + V'seABCD + i^i^BCD (39) 

with 

V^o = M: V^i = 2\/277i, -02 = 6^, 
•03 = -2\/2^, -04 = JI. 



In order to show the equivalence of the expressions psp and p9p one makes use of the spinorial 
identities given in appendix [Xj Following the notation of section 12.31 one writes 

■4'ABCD = ^iBCD + i'ABCD + i^ABCD + ^ABCD- 

The non-vanishing components of ^^abcd ^^^ given by 



^2 - ^3-^0 0- 



Those of i^ABCD ^^^ 



1 P / n , . r-i \ m D * ^ ,-1 rT-\ 1 



v-r = -,{P2 + iPi)r2°o - -yp^T2\ - -,{P2 - iPi)T2% 



4^ = --,{P2 iPi)TA ^-^PsT,\ + A(P2 + iPi)TA- 
P P P 

While in the case of ^^abcd '^^^ ^^^^ 

^3 = -^iQ2 - iQl)T2\ - ^Q3T2\ + 4(^2 + iQl)T2°2. 

p^ p^ p^ 

And finally for i^ABCD^ 

i'i - 4^(-Ji + iJ2)T2°o + ^-^J^T.'o - ^{Ji + -^J2)TA, 

^i = 0, 
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In the above expressions A, Pi, P2, P3, Ji, J2, J3, Qi, Q2, Qs G K- With regards to the part of 
ipABCD derived from the complex function one has that 

M2 = 2p9p(p9pA(^)) + X+X_A(^) - 2A(^) - pdpX^'\ 

with 

A = Xl'X = XlRcCX) + XllmCX) = A^^) + A^^). 

By direct inspection of the previous expressions of the components of ipABCD one has the 
following lemma. 

Lemma 4. The lifts to Ca of the solutions to the flat momentum constraint satisfy 

iPabcd[A] = p^^Eabcd[A], type{EABCD[A]) = p, 
'ipABCDiJ] = p^^Eabcd[J], typeiEABCoiJ]) =p+l, 
IpABCD [Q] = p^^Eabcd [Q] , type(5ABCD [Q] ) = P, 
i'ABCDiX'^^A = p"^5asci3[A('^], type(SABCD[A('^]) = P, 
i^ABCoiX^^^/r] = p-~^Eabcd[X^^^ /r], type{EABCD[X^^^ /r]) =p+l. 



Note that because iPabcd is totally symmetric, then its projection to Ba is associated to 

an /i-trace free symmetric tensor. In other words, borrowing the terminology of section 12. 4| 
c 

b 



<^ab = (Ja^Cr^^i^ABCD 



4.2 The non-conformally flat case 

A detailed discussion of the solutions, u*, of equation IT^ has been given in appendix [Bl From 
this analysis one can readily deduce the expansion types of the lift to Ca of the various parts of 
v\ 

Consider the vector u*[Re(A(^^)/r] produced by the seed V^y [Re(A'^^^)/r]. In appendix [Bl it is 
shown that under the assumptions of section 12.41 one has 



w'[Re(A(2))/r] = vi[Re(A('))/r] + w^[Re(A('))/r], 

with w*JRe(A(2))/r] e Qoo{Ba), wi[Re(A(2))/r] = 0{r^), and w^[Re(A(2))/r] e C°°{Ba). In what 
follows, the affixed [Re(A*^^-')/r] will be suppressed for ease of reading. Noting that because of 
Xiv"^ — r^u, u e C°°{Ba), one can write 

v\ = run^ + wvn} + Wfrf , 
with w a C°° complex function. Thus, one has that 

u* = rn^ + (wm^ + wm''). 
Accordingly, the lift of v"^ to Ca is readily found to be 

VAB = PUXAB + {WZAB + WTJAb) + v\bi 

with V, u, v\g of expansion type p + 1. Write uab — uxab, wab = wzab and wab = wjjab- A 
calculation renders 

D(ABVcD) = UX(abXcD) + PD(ABUCD) + D(abWcD) + D^abWcd) + D^abVcd) 

= P^^ ( P^ux^abXcd) + P^Di^abUcd) + P^D(^abWcd) + P^D(^abWcd) + P^D(^abVcd) j • 
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From the discussion in the previous paragraph it follows that type(p^M) = p—2, type{p^D(^jiBUcD}) 
p — 2, type{p^ D(^abWcd)) = P ^ 1 and type(p^Z?(^s''^c£>i) ^ -P ^ 1- Accordingly, one can write 

Using similar arguments one obtains the following 



Theorem 2. Under the assumptions of section ] 2. 4\ on the initial data, the lifts vab of the vectors 
v^ solving equation il^) are of the form: 



D(ABVCD}[A] = P~^<^\BCD[Jh iyP<^ ABCd[J]) = P, 

D^abVcdM = P'^'^abcdIQI iypeC^ABCDiQ])^?, 

D^ABVCD)[M^^'^)] = p-'nBCDiM^^'^)], type{nBCD[M^^'^)])=P-h 
D^ABVCD)[IH^^'^)] = P-'^"abcd[M^^'^)], type(ci>:^5c^[Im(A(i))]) =p- 1, 
DiABVcD)[M^^^^)/r] = p-^nBCD[M^^''')/r], type(<i>lsci?[Re(A(2))/r]) -p- 1, 
D^ABVcD)[lni{X^'^)/r] = p-3$^5^^[Im(A(2))/r], type{'^lscD[M^^^^)/r]) ^ P- 

5 Structure of the Weyl tensor on Ca 

The rescaled Weyl tensor plays a fundamental role in the discussion of the asymptotic properties 
of the gravitational field. Anticipating the discussion of section[6]an analysis of its structure near 
infinity at the level of initial data using the manifold Ca is now given. The discussion of this 
section is a generalisation of the analysis of section 4 in reference [H] . 

The spinorial counterpart of the rescaled tensor W^^Ap is given in terms of 

WaA'BB'CC'DD' — 4>ABCD^A'B'^C'D' + 4'A'B'C'D'^AB^CD, 

where (pABCD is the spinorial counterpart of a self-dual tensor. At the level of initial data, (J^abcd, 
is fully described in terms of the spinors wabcd and w'^g^jj — the spinorial counterparts of the 
spatial tensors Wab and w*^^ introduced in section [2. 1.11 One has that 

(j^ABCD — WABCD + iw^ABC£>- 

Using the conformal constraint equations (|8a p -(|8h p of section [^.1.11 and transcribing into space 
spinorial language one finds 

WABCD = ^^^D(abDcd)^ + ^~^SABCD + ^"^ ix^^ efX{abcd) - X^^ (abXcd)ef] , (43a) 

w*ABCD = -~''^^^^^D^(aXbcd)f, (43b) 

where xabcd is the spinorial counterpart of the second fundamental form Xij ■ For the class of 
initial data under consideration (maximal) X^^ ef = 0. If "ipABCD is the spinorial counterpart of 
the rescaled second fundamental form V'ij then 

XABCD = ^ IpABCD- 

As the magnetic part, w^^^jj^, is linear in ipABCD, it will prove useful to consider the following 
splitting oiw\g(.^: 



w. 



'Abcd[A. J. Q] + w*abcd[^'^''^ + w*abcd[>^'^''^] 



^ABCD 

In this last formula w\gij jj[\'^^^] denotes the part of w^^^^, obtained from il'ijiyS^^] while 
w\g(jgi[\''^''] is the part calculated from V'ij[A'"'^^]. 
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5.1 The massive and massless parts of (pABCD 
Using formula psp . the conformal factor $1 can be written as 

It-|2 



IJJ+\x\W)^' 
Let now 

In terms of the latter quantity one can calculate the massless part of 4>abcd- Namely, 

^ABCD — ^ABCD + '^ABCD' 



(44) 



with 



t>ABCD - ^ (u^D^abDcd)\x\^ - 4UD^ab\x\^Dcd)U ~ 2\x\^D^abDcd)U 

+6\x\^D^abUDcd)U+ \x\^U^sabcd] , 
\x\^ 2\/2 V2|a;P 

^ABCD = ~JJW^^^{AB^CD)EF + -JJ^D^ (^a\^\'^'^BCD)F - 4 ^3 D" ^^U i> bC D) F 
+ —[p-D (A^BCD)F- 



From here it is direct to obtain expressions for w'j^bcd ^^^ ''^abcd- One finds that 

I a; I lEF 
JJWV {ab 



w'abcd = 'f'ABCD ~ -rm'^'^''(AB'^CD)EF, (45a) 



w 



V ,1-^1 j.EF 



ABCD — VaBCD ' jjQ 



+ jk^ iABi^CD)EF- (45b) 



The massive part is given by 



^ABCD — ^ABCD + '^ABCD^ 



(46) 



^ABCD^r^i ~ t;^UWD^ab\x\^Dcd)\x\' + UW\x\D^abDcd)\x\' 



with 

1 / 3 

+2\x\{WD(^AB\xfDcD)U - 3UD(^AB\xfDcD)W) 

+2\x\^{- UD^abDcd)W - WD(^abDcd)U + 6D(^abUDcd)W + UWsabcd) 

+ \xW - 2WD(abDcd)W + QD(abWDcd)W + W^sabcd) 
^ ^'^ 

-ABCD U^{U + \x\Wf 

)EF 



■{W^''\x\W + lW'^\x\^W'^ + 2QU^\x\^W^ + lW'^\x\^W* 

+ m\x\^W'' + \x\^W')i;^^^j,Bi:cD)l 



^"^ {2U\x\W + \x\'W^)DP^j,\x\'^j,BCD)F 



U'^{U+\x\W) 



1 



+ L/3(^+|;'|V)3 (3^l^l^^^ + ^'l-l^ + l-l^^^)^Vt^'/' 



BCD)F 



AV2\x\'W ^f , ,^ , 4\/2| 



+ (£/ + NT4-)3 ^ (^l"l^^^^)^+(t7TNW (^^^^^^)^ 



V2| 



2 



C/2(L/+|a;|Ty) 
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{2U\x\W + WW^)D^^^j,4,BCD)F 



respectively, the time symmetric and non-time symmetric parts of the massive part of the Weyl 
spinor. 

It wiU turn out necessary to refine further the decomposition of the massless magnetic part of 
4>ABCD — which will be denoted by w^g^jj. Due to the linearity of the momentum constraint, 
it is possible to consider individually the various parameters in the second fundamental form. 
These will be denoted by w'^gQj) [A] , w*lg(jjj [J] , etc. 

From expression (|46|) for the massive part of the Weyl spinor one can see that unless the ADM 
mass of the data vanishes one has 'P'abcd ~ '^(l^^l"'^)- Thus, in order to discuss the behaviour 
of i^'abcd usar i one needs to introduce a suitable rescaling[^. To this end let n = \x\k' with 
K'{i) = 1 smooth. Consider the lifts to Ca of the spinorial fields 

4>ABCD = K 4>ABCD, 4>'aBCD = '^ 4>'aBCDi <PaBCD = '^ 'PaBCDj 

and so on. Let (J)abcd denote any of the aforementioned spinorial fields. From general principles 
one would expect 4>abcd to be of expansion type p + 2. That is, its essential components (f>j, 
J = 0, . . . , 4, have normal expansions near 2^ of the form 



J2 Yl J2-^^l^P'^1^kT^q''q+3-2P''^ 
p=0q=\2-j\ k=0 ^' 



with <f>j,p-2q,k G C. The symbol ^ is to be understood in the sense described in section 13.31 It 
turns out that the normal expansions have a more restricted form. The following generalisation 
of parts (i) and (ii) of theorem 4.1 in [11 will be proved. 

Theorem 3. The analytic lifts (J)abcd, 4''abcd '^'^'^ 'i'ABCD ^'^ ^a have expansion typep, whereas 
w'^g(jjy[A, J,Q] is of expansion typep— 1. In addition: 

(i) The expansion coefficients 0' 2q k °f ^'j ~ 'f'^ABCD) ■ satisfy 

^0,p;2p,fc = ^4,p;2p,fc' P = 0,1,2,..., fc = 0, . . . , 2p. 

(ii) The expansion coefficients w'- jr, k satisfy the antisymmetry condition 

^0,p;2p,fc — ~^4,p;2p,fcJ P = 0, 1, 2, . . . , k — 0, . . . , 2p. 

(Hi) The expansion coefficients w*' 2pfc[^jQiR-e(A)] satisfy the antisymmetry condition 

Wo[p-2pA^, Q, Re(A)] = -wl[p,2pA^^ Q' Re(A)], p - 0, 1, 2, ... , k = 0, . . . ,2p. 
(iv) The expansion coefficients ui*' 2r, ki'^i^'^'^i^)] satisfy the symmetry condition 

'*o!p;2p,fc[^)Im(A)] = W4'p.2p,fc[J,Im(A)], p = 0,1,2,..., k = 0,...,2p. 

The proof of the various parts of the theorem will be given in the following subsections. It 
consists, essentially, of an analysis of the various terms in expressions (|44| and (|46p of the massless 
and massive parts of the Weyl spinor. 

5.1.1 Proof of the part (i) of theorem[3] 

Recall the split 4>\gQ]j = Vabcd + 'tabcd introduced in equation (|46p . Consider first the term 
TABCD = '^^'/'abcd- This term coincides, formally, with the time symmetric 4>'abcd discussed 
in reference [11]. Note however, that in that refrence, W is the lift of an analytic function on 
Ba, while in the case treated here it is the lift of a function belonging to E°°(Ba)- However, by 
virtue of lemma[2]one has that type(iy) = p, and hence the argument in [11], which only requires 



^It is worth noticing that on the other hand, under suitable assumptions 4>'abcd ^'^ ^'^ analytic spinor on Ba 
These considerations will be retaken in section |9] 
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W having this expansion type (and not the analyticity) can be reproduced. This will not be 
repeated here. One obtains 

type(fec£,) =p, 

and the symmetry 

'^0*p;2p.fc = <^4*p;2p,fc. P= 0,1,2,..., fc = 0, . . . , 2p. 

Now, consider the term 4>Xbcd- From lemma|3]and theorem[2]it follows that the lift, 4'abcd 
to Ca of Tp^j = V'u [^j Q, J, A] + [Cv)ij satisfies ^Pabcd = P^^^abcd, with ^abcd a spinorial field 
of expansion type p+ 1. Hence, it follows that 4'^^ ij^gi^cD^EF = P^^'^^^ iab'^CD)ef^ where the 
spinor '^^^ iji^g^cD)EF has expansion type p + 2. Recall that multipHcation by a scalar function 
of the form \ + 0{p) does not change the expansion type or symmetries of the various terms. 

Using the observations in the previous paragraph one has that 

^^^^^ {W'pW+lW^p^W^+20U^p^W^+lW^p^W^+Wp^W''+p^W^)i^''^^j,s4cD)EF, 

— which is essentially the lift of \x\^ times the terms in the fifth and sixth lines of formula (|46)) — 
has at most expansion type p — 2. Similarly, 



lifts to Ca as 

which can be checked to have expansion type of at most p — 1. Now, 

lifts to 

hence, using that in the cn-gauge one has that type(C/ — 1) = p — 1 one concludes that the whole 
term has at most expansion type p ~ 2. The term 



\x 



5 



{u + \x\wY 

lifts to 



D ^^A\A'4^BCD)F, 



Y^X (^A^BCD)F, 



{u + \x\wy 

which can be seen to have at most expansion type p — I. The lift of 

is given by 

Now, DabW has expansion type p + 1, and hence D^,^W^bcd)f has expansion type p + 2. 
Accordingly, the overall expansion type of the term is p — 1. Finally, the lift of 

^^jJ^^(2C/|x|Ty + \xfw')D^^^i,BCD)F 
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is given by 
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Hence, noticing that D^ ,^^bcd)f f^nd x^ iji^^BCD)F have both expansion type p + 1 and p + 1 
one concludes that the overall expansion type of the term is at most p — 2. Thus, one has that 

type(0^BC£,) =P-1, 
and hence the symmetry 

^0,p:2p,k = ^4,p;2p,k^ P= 0,1,2,..., fc = 0, . . . , 2p, 

holds trivially. This proves part (i) of theorem [31 

5.1.2 Proof of the part (ii) of theorem [3] 

Recall now the expression (|45ap for w^scd ■ '^^'^ ^^^^ ^'^ ^o. of (p^bcd i^ identical to that of the 
time symmetric 4>'abcd discussed in [llj . The argument used in that reference to analyse the 
expansion type and symmetries of the time symmetric (J^'abcd uses in an essential manner the 
analyticity of the conformal metric hij to introduce a complex null cone formalism. To this end a 
3-dimensional complex analytic metric manifold (B, h) was introduced. In this setting h defines a 
complex valued non-degenerate scalar product. The complex manifold {B,h) contains {Ba,h) as 
a real Riemannian subspace. The spinor-dyad bundle SU{Ba) has a complex analytic extension 
to a bundle SL{B) of spin frames on B with structure group SL{2, C). In the same way that the 
metric h is analytically extended to B, one can also analytically extend the solder and connection 
forms. Crucial is now to consider the complex null cone. A/", generated by the geodesies through i. 
As one is restricted to work with analytic functions, one considers the analytic conformal factor 
il' obtained setting VF = in formula p^ . Because il' = p^/U, one then has that the null cone 
Af corresponds to the locus of points on B such that fl' . A suitable frame (i.e. coordinates and a 
frame) adapted to the geometry of A/" can be introduced on B. This results in a formalism which 
allows to calculate at the point i. A similar formalism has been used with great effect to discuss 
the convergence of multipole expansions of static spacetimes — see [13j . The complex null cone 
formalism is a powerful machinery to calculate the properties of the expansions of 4>'abcd ■ ^^ill 
the required calculations extend over 10 pages; it will be omitted. The key results are one obtains 

type(^ABC£)) =-P> 
and 

'f>Q,p;2p.k = -4>Xp;2p,k, P = 0, 1, 2, . . . fc = 0, . . . , 2p. 

Remark. This is the only part of the proof of theorem [3] where the analyticity of hij is used 
in an essential way. Still, it is expected that a similar result would follow in the smooth setting. 
This however, would involve lengthy induction arguments which will not be considered in this 
article. 

To conclude the analysis of w'j^bcd consider the lift to Ca of 



'^ F lEF I 

—(j&-i' (ABi^CD)EF- 



The lift is given by 



'''^^!$si=^,,„$. 



~fj&~^ (AB^CD)EF- 

As already discussed, ^abcd has expansion type p+ 1, and hence ^^^ /ab^cd)ef has expansion 
type p -\-2. Accordingly, the expansion type of the whole term is p — 1. 

Hence, the part of w', 2p k coming from ijjabcd satisfies trivially the required antisymmetry 
condition. This concludes the proof of part (ii) of theorem [3l 
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5.1.3 Proof of part (iii) and (iv) of theorem [3] 

Notice that using equation (|45bp one has 

WABCD = JpD (^|x| 1pBCD)F ~ ^ JJS ^ (aU'<PbCD)F + ^2 ^ {aVbCD)F- 

As in the proofs of parts (i) and (ii) consider, one by one, the hfts of the various terms on 
'^*ABCD ~ '^^'^ABCD- ^^*^ notes that the lift to Ca of 

-YpD^(A\xf^BCD)F 

is given by 

^X^iA'^BCD)F, (47) 

which has expansion type p as ^abcd has expansion type p + 1. Multiplication by the lift of k'^ 
will not alter the expansion type. The lift to Ca of 

Jp-D^(AU'4'BCD)F 

is given by 

— D (jJJ(^BCD)F- 

Now, using lenima[T]one can conclude that in the cn-gauge ty]ie{D abU) = p. Hence 

type(i:'^(^C/$BC£,)^) ^p+l. 
Accordingly, the overall expansion type of the term is p — 1. Finally, the lift of 

JpD^{Ai'BCD)F 

is given by 

2 

- 3^2;^(A<i'BCD)F + ^i?V*scD)F. (48) 

The expansion type of the first term of the latter expression has already been discussed — 
it is, modulo a constant, the same as the term (H7|) . To analyse the second term, note that 
although Dab^cdef has expansion type p + 2 if ^cdef has expansion type p+ 1, one finds that 
^^(A^BCD)F has expansion type p + 1 — for this, expand Dab^cdef in terms of symmetric 
irreducible terms and then contract indices. Multiplication by p^ renders a term with expansion 
typep- 1. 

Summarising, it has been found that the only terms in w'^bcd contributing towards an expan- 
sion of type p are proportional to pU^^x^ ,j^^bcd)f- The spinor ^abcd is totally symmetric. 
Hence one can write it as 

^ABCD = '^Q^°ABCD + ^1^\bCD + ^'^^'aBCD + ^se^BCD + ^^^ABCD- 

Furthermore, noting that the following relations hold 





r^F _ ^ T-F A - A 

^ (A^BCD)F — p:'^ABCD^ ^ (A'-BCD)F ~ „ /^<^ABC£li 


(49a) 




■^ (A'^BCD)F — "' 


(49b) 




F3 3 FA 4 
^ {A'^BCD)F — „ fr^ABCD^ ^ {A'^BCD)F /K^'ABCDy 


(49c) 


one finds that 






^F ^ 


1 ^_,0 , 1 ^_1 1 ^_,3 1 ^.,4 
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Now, ^ABCD has expansion type p + I, that is 

oo p+1 ^"^ 1 

p=0q=max{|l-j"|,2-p} fc=0 ^ 

A quick comparison with formulae (|42ap - (|42ap shows that the expansion coefficients of $o and 
$4 are essentially those of fi2 and /l2- Now, from lemma |4] and theorem [2] it follows that all the 
contributions of type p in ^abcd come from iJABColJ, ^^^^r]. Recall that A^^^ — A^^^ where 
A"^ is a complex function with real and imaginary parts in C°°{Ba)- Accordingly, write 

oo p+1 2q 

^^'^ /P-T.T.J2 ^fp-M-M + 9P-M.M) T,, V2/' (50) 

p=l q=2 k=0 

where fp-2q,k, gp;2q,k G C. The coefficients fp:2q,k satisfy the conditions 

/p;2g,fe = (-l)'''^Vp;2g,2?-fc> P=l,---, 9 = 2, ...,p+l, k = 0,...p, 

so that they are associated with the real part of A*^^' while the coefficients gpaq.k satisfy the 
conditions 

ffp;2g,fc = (-l)''+'^+^gp;2g,2q-fc, P=l,---, 9 =2,..., p+1, /s = 0, ...p 

SO that they are associated with the imaginary part of A^^^ . A direct calculation reveals that 

oo p+l 2q 

^2q q-2l^ 

oo p+l 2q 



oo p+l 2q 

M2[Re(A(2))/p] ^Y^Y^Y.- {2p{p ~ l)P, - (g - l)(g - 2) - 2) /p;2,,fer2, 



p=l g=2 fe = ^ 



M2[Im(A(^))/p] - -EEE7Av5p;2,,.T,/,_,p^ 

p=l g=2fe = ^-^ '^^ 



with 



(iq = ^{q~l)q{q+l){q + 2). 
From here another short computation using the rules (|49a[) - (|49cP shows that 

$0,p;2p+2,fc[Rc(A(2))/p)] ^ _$4p.2p^2^^[Rc(A(2))/p], 
$0,p;2p+2,fe[Im(A(^))/p] = $4,p;2p+2,fc[Im(A(2))/p]. 

The latter relations imply the antisymmetry and symmetry conditions, 

wio,p;2p,fe[Re(A(2))/p] = -w;4^p.2p,fe[Re(A(2))/p], 
u'o,p;2p,fc[Im(A(2Vp)] = -W4,p;2p,fc[Im(A(2))/p]. 

Finally, a similar argument with pf/^^a;^.^V'BC£i)F shows that ^o[J] = ^4[J] so that the sym- 
metry conditions are trivially satisfied. This proves points (iii) and (iv) of theorem [3 □ 

5.2 Some further results 

Besides the information provided by theorem [3] about the properties of the expansions of (J)abcd 
and derived spinors, a couple of further results will be required. The following spinor will be used 
in the sequel 

CABCD = ''^V2D^{aXbcd)f- (51) 

The following lemma is a direct consequence of the proof of theorem [3l 
Lemma 5. For the class of data under consideration 

type{cABCD) =P+ 1. 
26 



In accordance with the previous result one has that the essential components Cj of cabcd are 
of the form 

oo 2p+2 2q 

Cj "" 2^ Z^ Z^ ~\'^0-VM-kT2q q-2+jP ' 

p=0g=max{|l-j|,2-p} fc=0 

Now, recall the expansion ([50)1 of A^^^/p. The coefficients fp-^2q,k associated to the real part of 
A*^^^ satisfy the following lemma. 

Lemma 6. For the class of initial data under consideration: 
(^ 

Kv;1vA^<^ I P)\ = ^ /p;2p+2,fc = 0. 

(n) 

Cj\p;2p+2,fc[^e(A(2V/0)] = ^ fp;2p+2,k = 0. 

Proof. The proof of (i) follows directly from the analysis carried in the proof of parts (iii) 
and (iv) of theorem [31 The proof of part (ii) follows directly from (i) by noticing that w'^^^jj = 
-Q~'^CABCD- n 

Finally, to state the main result of this article one needs the following theorem, which for the 
class of data under consideration, generalises part (iii) of theorem 4.1 in [TT] . 

Theorem 4. For the class of initial data under consideration one has that: 
fi) 



if and only if 

(n) 

if and only if 



i"o,p:2p,fc = 0, p = 2,3, ..., k = 0,...,2p, 

D(ApBp ■ ■ ■ DA^BtbABCD){i) =0, p = 0, 1, 2, . . . . 

<p;2p,fc[^e(A('Vp)] - 0, p = 2,3,..., fc = 0,...,2p 



D(A,B, ■ ■ ■ i?AilBlc^Bc)''[^e(A^'VO](*) = 0, p = 1, 2, 3, . . . . 
Remark 1. Note that the class of data under consideration automatically satisfies 

CABCD{i) = 0. 

Proof of theorem [4l From the proof of part (ii) of theorem [3] one has that 

type(0ABCD) = P' type('^Ai3CD) = P - 1- 

Hence, itn order to prove part (i) of theorem[31one has to consider only the expansion of the term 
TABCD which is formally identical to the spinor (J>'abcd considered in |11[ under the assumption 
of an analytic conformal metric. Thus, (i) follows from the analysis based on the complex null 
cone formalism given in the aforementioned reference and which will be omitted here. 
In order to prove part (ii) notice that consistent with lemma [5) 

^3 Z^ Tj\ ^3 P ' 

p=0 ^' 

with 

2p+2 p 2q 

c/ = ZJ '^^■P'2P+2.fc^2p+2 p-l+j + Z2 2_/'^J'P;2g,feT'2g 9+i' P=l,2,..., 

k=0 g=max{|l-j|,2-p} fe=0 
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where from general considerations about normal expansions at X — see section l375l and in particular 
the expression (|37p — one has that 



Cj,p:2p+2,k= Kpj^kD(^ApBp---DA,\E\Cj^SC), (*)' J = 0' ■ ■ ■ ' 4, fc = 0, . . . , 2p, 

and Kpj^k are some constants. Part (ii) now follows from lemma [6] □ 

Remark 2. It is perhaps worth noticing that cabcd being of expansion type p + 1, then 
Cj,p;2p+4,k — and accordingly, 

D(ApBp ■ ■ ■ DAiBiCABCD){i) = 0. 

6 The spacetime Friedrich gauge 

The final aim of the framework developed in pTj is to gain control over the evolution of the 
gravitational field in a neighbourhood of spacelike infinity which extends to null infinity. This 
problem is generically known as the initial value problem near spatial infinity. In slight contrast 
to the analysis of the non- linear stability of the Minkowski spacetime — see e.g. [2j |13], the aim 
of the initial value problem near spatial infinity as discussed in [llj is not only that of showing 
that the outgoing null geodesies starting close to spatial infinity are complete, but also to analyse 
under which conditions on the initial data, the resulting spacetime will admit a smooth conformal 
extension through null infinity — and hence giving rise to an asymptotically simple spacetime. 

In the standard representation of spatial infinity as a point, the direct formulation of an initial 
value problem for the conformal field equations with data prescribed in a neighbourhood, say Ba 
of infinity, renders a problem which although local is also singular — this can be seen in a very 
poignant way by considering the expressions for the massless and massive parts of the Weyl spinor 
discussed in section [5] 

The formulation of the initial value problem near spatial infinity presented in [llj employs 
gauge conditions based on timelike conformal geodesies. The conformal geodesies are autoparallel 
with respect to a Weyl connection — i.e. a torsion-free connection which is not necessarily the Levi- 
Civita connection of a metric. An analysis of Weyl connections in the context of the conformal 
field equations has been given in lOj. In terms of this gauge based on conformal geodesies — 
which shall be called the Friedrich gauge or F-gauge for short — the conformal factor of the 
spacetime can be determined explicitly in terms of the initial data for the Einstein vacuum 
equations. Hence, provided that the congruence of conformal geodesies and the fields describing 
the gravitational field extend in a regular manner to null infinity, one has complete control on 
the location of null infinity. In addition, the F-gauge renders a particularly simple representation 
of the propagation equations. Using this framework, the singular initial value problem at spatial 
infinity can be reformulated into another problem where null infinity is represented by an explicitly 
known hypersurface and where the data are regular at spacclike infinity. The construction of the 
bundle manifold Ca and the blowing up of the point i e Sa to the set l" C Ca briefiy described 
in section [3] is the first step in the construction of the regular setting. The next step in the 
construction is to introduce a rescaling in the frame bundle so that fields that are singular at T^ 
become regular. The required rescaling has already been hinted in theorem [3] where the regular 
spinorial field <J)abcd — k^4>abcd has been considered instead of the singular 4>abcd- 

6.1 The manifold Ma,K 

Following the discussion of [ll] assume that in the development of data prescribed on Ba the 

timelike spinor r introduced in section [3] is tangent to a congruence of timelike conformal 

geodesies which are orthogonal to Ba- The canonical factor rendered by the consideration of this 

congruence of conformal geodesies is given in terms of an affine parameter r of the conformal 

geodesies by 

/ K^r^N 20. 

e = K-^fl 1 —], with uj^^ (52) 
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where fl = i?"^ and -d solves the Licnerowicz equation (|13p . The function k > — which will be 
taken to be of the form k — k' p, with k' smooth, k' {i) = 1 — expresses the remaining conformal 
freedom in the construction. Consistent with the scalings 6a ^^ i^^^'^^a induced by the function 
At one considers the set Ca,K — K^^'^Ca of scaled spinor dyads. Furthermore, define the bundle 
manifold 



K{q) K{q) 

which, assuming that the congruence of null geodesies and the relevant fields extend adequately, 
can be identified with the development of Ba up to null infinity — that is, the region of spacetime 
near null and spatial infinity. In addition to Ca,K, one defines the sets: 

I={{T,q)^Ma,.\p{q)=0, |t|<1}, 
I± = {(T,q) e Ma,. I p{q) = 0, r - ±l}. 



K{q) 

which will be referred to as, respectively, the cylinder at spatial infinity, the critical sets (where 
null infinity touches spatial infinity) and future and past null infinity. In order to coordinatise the 
hypersurfaces of constant parameter t, one extends the coordinates (p, t^g) off Ca,K by requiring 
them to be constant along the conformal geodesies — i.e. one has a system of conformal Gaussian 
coordinates. 

6.2 The conformal propagation equations 

On the manifold Ada.n it is possible to introduce a calculus based on the derivatives dr and dp 
and on the operators X+, X^ and X. The operators dp, X+, X^ and X originally defined on 
Ca can be suitably extended to the rest of the manifold in a standard way. A frame caa' and 
the associated spin connection coefficients Taa'bc of the Wcyl connection V will be used. The 
gravitational field is, in addition, described by the spinorial counterparts of the Ricci tensor of 
the Weyl connection, Qaa'bb', and of the rescaled Weyl tensor, (jiAscD- In order to describe the 
conformal propagation equations consider the vector 

^ — {caBi^ABCDiQaBCd)i (t> — {(pABCo), 

where cab, ^abcd, Qabcd are the space-spinor versions of the spacetime spinors caa' , ^aa'BC, 
^AA'BB' ■ The relation between space-spinors and spacetime spinors is implemented by means of 
suitable contractions with the spinor t^^ . This will not be elaborated further — the interested 
reader is referred to [El [6l [11]. The explicit form of the propagation equations will not be 
required. Only general properties will be used. Schematically one has 

drV ^ Kv + Q{v,v)+L(j), (54a) 

V2Edr^ + A^''c'XBdp(l)^BiTABCD)(l), (54b) 

where E denotes the (5 x 5) unit matrix, A^^c'^^ are (5 x 5) matrices depending on the coor- 
dinates, and B{Tabcd) is a linear (5 x 5) matrix valued function with constant entries of the 
connection coefficients Tabcd- In addition to the above propagation equations it is essential to 
consider the following constraint equations derived from the Bianchi identities: 

F^^'CXsOpCJ) - HiTABCo), 

where now F^^c'^^ denote (3 x 5) matrices, and H(rABCD) is a (3 x 5) matrix valued function 
of the connection with constant entries. 

Equations (|54ap and (|54b|) can be casted as a symmetric hyperbolic system on a neighbourhood 
Af C Ma,K. of the initial hypcrsurface. Hence, given data that extends smoothly to l", one obtains 
a unique smooth solution on J\f. The metric can be seen to degenerate as p -^^ 0. A consequence of 
this is the fact that A^^ c\g = onXnA/" — i.e. the coefficient associated with the p derivatives 
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in the Bianchi propagation equations (j54bp . The restriction of the propagation equations to X 
imphes an interior system on X which determines v and (j) oti X uniquely from the restriction 
of the initial data to X" . Differentiating the propagation equations repeatedly with respect to p 
and restricting the result to X one obtains a hierarchy of interior equations for v'^^^ — d^v\-i from 
where it is possible to determine a formal expansions 

of the solution at X. The calculation of the lowest order terms u^"-' shows that the matrix A-^^c^g 
is positive definite on X and extends smoothly to X , but it looses rank there. The entries of 
the vectors v^^' and (p^^' on X can be seen to have definite spin-weights and hence admit very 
particular expansions in terms of the functions T- ^ — as described in section 13.51 The use of 
these time dependent normal expansions reformulates the problem of calculating the vector v^^' 
into a problem of linear ordinary differential equations. Once the entries of the vectors v^^' and 
(f>^P' have been expanded using the functions T- j, the interior equations are reduced to systems 
of ordinary differential equations for the r-dependent expansion coefficients. The task of solving 
these equations reduces, in the end, to solving a hierarchy of ordinary differential equations of 
the form 

Va^CaVa+ba, (55) 

where Cq, is a 2 x 2 matrix and j/^ and ba are 2x1 column vectors. The components of j/q consist 
of pth-order p derivatives of certain components of the Weyl spinor <f>ABCD- In equation (|55|) note 
the presence of a multi-index a = (p, q, k) indicating the order, pi^, in the expansions and to which 
harmonic T2 ; the term is associated. In what follows, in order to ease the discussion, the multi- 
index a will be, sometimes, suppressed. For a given multi- index a = {p,q,k), the components 
of the vector b are calculated from the lower order solutions w^"^-* and 4''"'^\ < q < p — 1. The 
solutions to these equations can be written in the form 

yiT)=X{T)X-\0)yo+X{T) f X-\s)b{s)ds, (56) 



where yo — y(0) and X{t) — again suppressing the relevant multi- index — denotes the fundamen- 
tal matrix of the system of ordinary differential equations. The matrices X have been explicitly 
calculated in [11] . For increasing p the explicit expressions for b become more complicated. One 
can, nevertheless, implement the aforediscussed procedure in a computer algebra system — see 
e.g. [mEl]. Now, due to the degeneracy of A^^c\g on J, the ordinary differential equations 
for the vector y are singular at r = ±1. As a consequence, for a given p, there are certain choices 
of the multi- index a, = {p,p,k), k — 0, . . . , 2p, of for which the fundamental system develops 
logarithmic singularities at r — ±1. For all other allowed values of the multi- index, the funda- 
mental matrices can be written explicitly in terms of Jacobi polynomials. It can be seen that the 
vector b for the corresponding multi-index a* vanishes. Hence the logarithmic divergences cannot 
be cancelled out by the integral in ([56|) . Thus, in general the functions -u'^^ and 0*^^^ develop 
logarithmic singularities for all p > 2 on X. 

6.3 Regularity conditions 

Given the aforementioned state of affairs, can one find conditions on the initial data so that the 
functions j/q extend smoothly to X^l An inspection of the term X{T)X~^{0)yo for the singular 
multi-index values a* shows that there are conditions on the initial data for which the observed 
singularities at X do not arise. This analysis is completely general as long as the initial data is 
expandable in powers of p near X". 

With the aim of formulating the main result of this article, it will be necessary to provide 
a more precise description of the conditions 00 the data mentioned in the previous paragraph. 
Following the conventions of section [5] consider the spinor (pABCD = k'^'Pabcd, and let 
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Observing theorem [3] one expands 

p 2q 

^f = 2^ 2^ %:P;2g,fc^2q q-2+j' 
q=\2-j\ k=0 

with r-dependent complex coefficients aj^p-2q,k- The coefficients aj,p-^2q.k can be determined along 
I from their value at 1° by solving a set of transport equations implied by the conformal field 
equations as described in section [6] The analysis of [Tl] implies the following result. 

Theorem 5. For the class of data under consideration one has that 

ao,p;2pAr) = (1 - rr'il + rf-' (^o,. + C,^k [ (, ^ ^),,f(, _ ^^,^3 
a4,;2p,.(r) = (1 + rY^\l rf-^ (co,. + Ci,. ^^ ^, _ .y-'l, ^ ,y+, ) - 

wjt/i Co,fc and Ci^k constants. In particular, aQ^p:2p,k and a4^p-2p,k extend analytically through 
T — ±1 if and only if 

a0,j9;2p,fc(0) = a4,p;2]9,fe(0), 

with k = 0, . . . , 2p. 

7 The main result 

In [11] it was shown that if one restricts the attention to the class of time symmetric data 
with a conformal metric that is analytic in a neighbourhood of infinity these conditions can be 
reformulated in terms of the vanishing of the Coton tensor and its symmetrised higher order 
derivatives at infinity. This condition happens to be a purely asymptotic condition on the freely 
specifiable data. 

The main result in this article is a generalisation of the result described in the previous 
paragraphs to the class of data with a non- vanishing second fundamental form. This main result 
brings together the discussion in section [5] (on the structure of initial data for the Weyl tensor 
near infinity), in particular of theorems [3] and |31 with the discussion of section [S] (on the regular 
finite initial value problem at spatial infinity) leading to theorem [5] 

From theorem |3] it follows readily that 

a0,p;2p,fc(0) = 4>^p.2p,k + ^ip;2p,k + Wo[p;2pA^'^i>'''^^ / P)] + W^ap;2p,fc [Im(^'^Vp)] > 

a4,p;2p,fc(0) = (i>^p;2p,k - w'o,p;2p,k " w^S'p;2p,fc [R.e(A(2Vp)] + ^o'p;2p,fc[Im(A(2Vp)]- 
Hence the condition 

a0,p;2p,fe(0) = a4,p;2p,fc(0), fc = 0, . . . , 2p, 

of theorem[5l using that wabcd is associated to a real spatial tensor whereas w\gujj is associated 
to an imaginary one, implies 

4,p;2p,fc = 0, wl[^,2p.k [Re(A(2) /p)] = 0, 

on Ca independently of the choice of k. These last two conditions can be readily reformulated in 
terms of the tensors bABCD and cabcd using theorem |4l In this way one obtains the main result 
of the article. 

Theorem 6. For the class of data under consideration, the solution to the regular finite initial 
value problem at spatial infinity is smooth through T^ only if the conditions 



D{ApBp---DAiBibABCD){-i) ^0, P = 0, 1, 2, . . . (57a) 

D^A,B,---DA,\E\c^sc)''^[X^''^]{i) = 0, g = 0,1,2,... (57b) 



are satisfied by the free initial data. If the above conditions are violated at some order p or q, 
then the solution will develop logarithmic singularities at I . 
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Remark 1. This result is a non-time symmetric generalisation of theorem 8.2 in [llj . 

Remark 2. The result in theorem [6] which is expressed in the language of space-spinors can 
be readily reformulated in terms of spatial tensors. Using standard transcription rules one finds 
the tensorial version of the main theorem given in the introductory section. 

8 Extensions 

In view of the main result [5] it is natural to ask whether the assertions are true for more general 
classes of initial data. In particular, it is important to see how stationary initial data fits into 
this picture. In [9] it has been shown that the condition (|57aP for the spinor Iiabcd is satisfied 
by static data — the condition (I57b|l is in this case satisfied trivially. 

In [3] it has been shown that there is a gauge for which in a neighbourhood Ba of infinity, the 
conformal metric of stationary data is of the form 

h^j^li^+r^h!^. (58) 

with hi-' and /i^- analytic. This type of conformal metrics is not smooth analytic at i — i.e. at 
r = {]. In fact, it follows directly that it is C^'". For the tensor ijjij it was shown that in that 
gauge stationary data are such that 

V'y-r-Sv.y^+r^Vlf, (59) 

with ■0J- ■ = 0{r'^) and ^j- — 0{r^) analytic. A choice of hij and Vij consistent with (|58|) 
and (|59(l provide a natural ground for generalising the results in this article. Unfortunately, this 
analysis is complicated by the fact that up to date there is no generalisation of the analysis in [5] 
for conformal metrics which are non-smooth. 

9 Conclusions 

In order to bring further context to the content of the main theorem[6l some connections with the 
construction of data for purely radiative spacetimes are raised: in reference [9] it has been shown 
how time symmetric Cauchy initial data with vanishing mass — so that O = fi' in the notation 
of section [5] — can be used to obtain analytic data at past null infinity. The crucial point in that 
analysis was to obtain conditions on the Cauchy data so that the spinor (Jjabcd = wabcd is 
analytic at i. The Cauchy-Kowalevskaya theorem can then be used to obtain analytic data on 
past null infinity close to i~ (purely radiative data). The conditions are given by 

D(A,B, ■ ■ ■ DAiBibABCD){i) =0, ? = 0, 1, 2, . . . , 

that is, the same as condition (|57ap in theorem [Bj The way the analysis in reference [S] can 
be generalised to the case of non-time symmetric Cauchy data has been briefly discussed in 
the concluding remarks of [H]. That analysis will not be repeated here, but under the extra 
assumption that the Cauchy data has no linear momentum — so that Xij = ^('') — j necessary 
and sufficient conditions for (J)abcd = wabcd + i^ABCD ^"^ ^^ analytic at i are in addition to 
(I57a|l that 

DiA,B,---DA,B^)CABCDil)^Q. p = 0, 1, 2, . . . . (60) 

Note that condition (I60|l is much stronger than condition (j57bp . Indeed, from the discussion in 
subsection 13.31 one has that 

^(ApBp • ■ • Da^Bi)CABCD = ^ D(^ApBp ' ' ' Dai\E\Cabc) = 0, 

but not conversely. Furthermore, condition (|60p involves the whole of ipij and not only the part 
depending on Rc(A^^Vp)- 
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A Spinorial identities 

The following identities will be used throughout the main text: 

_ r, 2 _ 1 4 

X(ABXcD) — '^^ABCDi y(AByCD) — -^^ABCD^ 

_ 3 _ 1 2 

X(AByCD) — -^ABCDi y(ABZCD) — -^^ABCD^ 



. . _lo 



Also, 



and 



^(ABZcD) — ^ABCDi Z(ABZcD) — T^^ABCD- 



Xa^XbF = -^^AB x/'ubf = JkVAB, x/" ZbF = -j=ZAB, 



F _ ^ F 1 _ 1 1 

^ {A^BCD)F — "Tk^ABCD^ ^ {A^BCD)F — ~Jk^ABCD^ 

X {A^BCD)F — ^1 

F 'A 3 FA 4 

^ {A^BCD)F = Ik^ABCDi X (^a^bCD)F = Ik^ABCD- 

B Technical results concerning solutions of the constraints 

A detailed analysis of solutions to the vacuum Einstein constraint equations such that near 
spacelike infinity they admit asymptotic expansions of the form 



hj r 


/ 2m> 
V ^ J 


Ik 
fe>2 


Xij '' 


fe>2 





where h^^ and x*^- are smooth functions on § has been given in reference [3- The results given 
in the latter reference will constitute fundamental building blocks of the analysis carried out in 
the present article. This appendix offers a summary of the results of reference [5] relevant for the 
analysis together with some extensions of their work as the class of initial data considered in the 
present article turns out to be more general than the one in reference [5]. 

B.l Solutions to the Hamiltonian constraint 

A function / e C°°{S) is said to be in E°°{Ba) if on Ba one can write / = /i + r/2 with 
/i, /2 S C°"{Ba), with r^ = SijX^x^ . One could define in an analogous fashion the spaces E^{Ba) 
and E'^{Ba)- Theorem 1 in [5] then states that: 

Theorem 7. Let hij be a smooth metric on S with positive Ricci scalar. Assume that ipij is 
smooth in S and satisfies on Ba 
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Then there exists on S ^ S \ {i\ a unique solution ^ to the Licnerowicz equation il3\) . which is 
positive, satisfies liiii,.^o r-d — 1 and has in Ba the form 

'&^-{ui+ru2), uieBa, U2eE°^{Ba), ui = l + 0{r^). 
r 

A class of symmetric trace-free tensors V'ij solving the momentum constraint and satisfying 
the condition r^ipijip^^ G E°°{Ba) has been discussed in section 4.3 of [S] and their asymptotic 
properties in corollary 5. Unfortunately this class of tensors ipij is not general enough for our 
purposes. Accordingly, further analysis is required. 

B.2 Solutions to the momentum constraint 

B.2.1 On the solutions to the flat space momentuni constraint 

As in the main text, denote by tpij a solution to the Euclidean momentum constraint and write 

i'ij = ipij [P] + ij'ij [A] + ipij [J] + ipij [Q] + -0ij [A] , 

where the [P] denotes the part of tpij depending on P*, etc. It can be seen that — see theorem 15 
in reference 0— that if rX e E°°{Ba) and F* = then r^-ip^j^P^^ e E°°{Ba). This result justifies 
the choice of A made in the main text — cfr. Ansatz (HI]). The /i-trace-free part 

of ipij will be used to specify the freely-specifiable data in the solutions to the momentum 
constraint. As mentioned in the main text, if r\ S E°°{Ba) then A = A*^^^ + yS'^^ /r, with 
^(1)^ ^(2) g C°°{Ba)- Accordingly, one can write without loss of generality 



A '^^ / ^i\ ■■■ik "^ ' * * *^ 7 ^ ^^ / Pii ■■■ik "^ 



k>2 k>2 

with aii-ifc , Pii---ik G C, symmetric. It is recalled that the above expressions are to be interpreted 
as 

m 771 

A«=^a.,....,x^^...x^'=-KAi,^\ A(2)=^A,....,x^^...x^'=-KAi,'\ 

fc>2 fe>2 

with A^^^ = air"") and A^^^ = o{r"') for all m. 

The expressions given in section 12.31 for the various possible solutions of the momentum con- 
straint are given in terms of its components with respect to a frame {n*, TO*,m'}. This choice is 
not ideal for the subsequent discussion of the solutions of equation (H^). Therefore, the expres- 
sions involving the parts of the solution arising from the complex function A are reformulated 
in terms of a Cartesian basis related to the normal coordinates {x*}. The lengthy calculations, 
which are omitted here, reveal that the tensor tpij [Re(A^^^)] is a sum of terms of the form 

.(1) o(l) ^(1) 



iXBy(Xi^...ij^ jX '''X X^Xj ^ „ r\,e(^G;j^ ...^^ jX •••X ^ij ' o '^(i^^\^j)ii---ip — i)'^ 



where Af.', bI,' and C^ are constants depending on fc > 2 such that the whole term is i5-trace 
free and (5-divergence free. The contributions due to the imaginary part of A^^^ render a tensor 
ipij[Ira{X^^')] which is a sum of terms of the form 

(1,0) n(i^2) 

kl T / 11 19 \ i\ ^, 



/ -L-'m kl T / \ i-i i™ 1 I ^-^rn kl -r / n 7 2 \ ii 



34 



where again D\. ' ' , D^, ' ,■■ ■ and Ej^ ' ,Ej^'' ,■■ ■ are constants making the whole term S- 
divergence free — note that the terms are 5-trace free by construction. Similarly, the tensor 
'0ij[Re(A*^^^)/r] is a sum of terms of the form 

^(2) ^(2) ^(2) 

-^Re{P^,...i,)x'' ■ --x'-x^Xj + -^Re{Pi,...^,)x'^ ■ --x'^S,, + -^a;(,Re(/3j),,...,^_Ja;*i • ■ ■ x^-^ , 

where A),\ B^, and Cj, are constants depending on fc > 2 such that the term is 5-trace free 
and 5-divergence free. In addition, one has that i/jij [Im(A^^-')/r] is a sum of terms of the form 

p(2,0) p(2,2) 



f^l Tm/'fl A^H . . . ^»,„-2 _L il!!!^ ^.^kl Ttth/'/? JiJ2^2;*i 



Xfee'='(Jm(/3,),,^...,^_Ja;^i • • -x^'-^ + ^:^a:fce'^-'(,Im(/3, 



where Z?^ ' , Dj^ ' , . . . and Ef. ,Ej. , . . . are constants making the term 5-divergence free 
— note that the individual terms are <5-trace free by construction. 
For future reference it will be convenient to write 

4[Re(A(i))] - r-5s,,[Re(A(i)], ^,,[Im(A(i))] = r-4s,,[Im(A(i))], 
^,,[Re(A(2))/r] = r-^E,,[Re{X^'^)/r], ^,,[Im(A(2))/r] = r-5s,,[Im(A(2))/r], 

with Sy[Re(A(i))], Sy[Im(A(i))], Sy[Re(A(2))/r], S,j[Im(A(2))/r] e C°°{Ba) —and actually 
polynomial. It is noticed that 

x*xJ5y[Re(A(i))] = r2s[Re(A(i))], a;'a:^S,,[Im(A(i))] = 0, 
a;*a;^Sy[Re(A(2))/r] = r2s[Re(A(2))/r], a;'a;JSy[Im(A(2))/r] = 0, 

for S[Re(A*^"'^-')], S[Re(A'^^)/r] some scalar functions. 

B.2.2 General issues concerning the solutions to the non-flat momentum constraint 



An analysis of the solutions to the elliptic equation for the vector u* — see equation p2)) — for free 
data of the form <I'ij[^, J, (3,Re(A*^^-'),Im(A^^'')/r], has been given in theorem 16 and of reference 

m- 

If the conformal metric hij admits no conformal Killing vectors on 5, then there exists a 
unique vector w* G W^''^, q > 1 solving equation P^ . If the conformal metric admits conformal 
Killing vectors one can guarantee the existence of a solution only if the constants A, ,P and Q* 
satisfy a particular relation. In order to discuss the asymptotic expansions of the vector u* recall 
the definition of Qoo{Ba) given in section [2.5.21 Furthermore, introduce for m e N, m > 1 the 
following real spaces of functions: 

Qra = {v' e C°°(M^R=^) I V' e Vra. v'x, = t" V with V G V-m-l) , 

where Vm denotes the space of homogeneous polynomials of degree m — that is, v G Vm if and 
only if w = Vi^.-.i^x^^ ■ ■ -a;'™, with Vi^...i^ G K totally symmetric. 

It turns out that a convenient way of grouping the different terms in the free specifiable data 
is the following: 

^y = i^y[A, J,Q] +^„[Re(A(i)),Im(A(2))/r] +^„[Re(A(2))/r,Im(A(i))]. 

We shall proceed to analyse the asymptotic expansions of the solutions to P^ implied by each 
of these terms. However, first an analysis of the source term D^<tij will be required. 
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B.2.3 Analysis of D*$y 

As seen in sections 12.41 and 13.61 , if the conformal metric hij satisfies the cn-gauge, then one can 
write 

hij = -S,j + r^h.j , h'^ = -% + r^h'^ . 

In addition, the technical assumption ([22|) is recalled: 6^^ hij — 0, Sijh^^ = 0. A calculation 
renders 

T/cy = {xihjk + Xjhik - Xk'hij) + -r^ {dihkj + djhik ~ dkhij) 

— -1- kij + r 1 kij • 

In addition, let 



One readily verifies that 

and that 

Now, recall that 

and that 
so that 



ij = n L kij, i ij = n I ki] ■ 



x'f\^ = r^h'^, x'f\^ = r^h\ xif\-^^ = r^hj, 



D'^,j = r^h'''^k^^J - h^'T\,4,i, - h^'T\^^a- (61) 

In the last equation it has been used that by construction 

5'''^k^^J = 0. 

Now, consider $ij[Re(A'^^^)/r], which is the most singular contribution of A to the seed tensor 
'hij. A direct calculation shows that r^h^'^dkipij is given by a sum of terms of the form 

j(2) o(2) (^(2) 

-^h'H,kReiP^,...,Jx'' ■ --x'-Xj + —^h'''S,,Re{P,,...,^_,k)x'' ■ ■ ■ x'-^ + -^h'''^5kiRe{(i^,...^,^_,J)x'' ■ ■ ■ x'-^ 
+ ^h'''5kjRe{P,,...,^_,,)x'' ■ • -x'"- + '^" ^'!~^^ Xj-fe"'Re(A,...,„_,fc,)a^'^ • --x'-'. 



Hence, if the technical condition ([22)1 holds, one can conclude that 

i?%[Re(A(2))/r] = l5,[Re(A(2))/r], 

with 

5,[Re(A(2))/r] e Qoo(S„), 5,[Re(A(2))/r] = 0(r^). 

A similar analysis can be carried out with the contributions due to Re(A'^^), Im(A'^^') and 
Im(A(^^)/7-. One obtains the following lemma. 
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Lemma 7. For the class of initial data under consideration expressed in the cn-gauge, and 
assuming that condition i22\) holds, one has that: 



1 



0\ 



D'<^,j[A] = -S,[A], S,[A] e Qoo{Ba). Sj[A] = 0{r , , 

D"^^o[J] = ^Sj[Jl Sj[J] e Qoo(6a), Sj[J] = Oir^), 

D'kM - ^S^m S,[Q] e QooiBa), S,[Q] = 0{r), 

i^^$,,[Re(A(i))] = i5,[Re(A(i')], ^,[Re(A(i))] e Q^Ba), ^,[Re(A(i))] = Oir'), 

D4,,[Im(A(i))] = l5,[Im(A(i))], ^,[Im(A(i))] G Qoo(S„), ^jMA^D)] - Oir^), 

i^4,,[Re(A(2))/r] = Ls,[ReiX^^^)/r], ^,[Re(A(2))/r] e Qoo(6,), ^,[Re(A(2))/r] - ©(r^), 

i^'^$,,[Im(A(2))/r] ^ l5,[Im(A(2))/r], 5,[Im(A(2))/r] e Qoo(Sa), 5,[Im(A(2))/r] = ©(r^). 

B.2.4 Free data V-ij [^, J, Q] 

The case has been discussed in corollary 5 of [5\ . Using the results of lemma [7] which hold for 
the cn-gauge and assuming condition (HH), a direct application of the techniques of section 4.2 of 
reference |S] render 

v^[A]=rvl[A]+v^,[Al 
vV]^lvl[J] + vl[J] 

v'[Q]^^vl[Q]+v',[Q], 

with 

vl[A]&Q^{Ba), v{[A] = Oil), vl[A]€C^iBa), 
v{[J]e QooiBa), v{[J] = Oir'), v'^[J]&C^iBa), 
vl[Q]e QooiBa), v{[Q]^Oir), v',[Q] e C°°iBa). 

B.2.5 Free data V'ii[Re(A(^)),Im(A(2))/r] 

This case has not been discussed in complete generality in [S] , but in view of the results of lemma 
[7] it is essentially a direct application of theorem 17 there. One obtains that 



w*[Re(A(i))] = rwi[Re(A(i))] +z;^[Re(A(i))], 
v'[Im(A(2))/r] = -t;i[Im(A(2))/r] +w^[Im(A(2))/r] 



with 



t;i[Re(A(i))] e QooiBa), v\[ReiX^-''^)] = Oir^), v',[ReiX^'^)] G C°°(S,), 
i;i[Im(A(2))/r] e goo(Sa), vl[ImiX^^'^ ) / r] = ©(r^), z;^[Im(A(2))/r] e C°°iBa) 

B.2.6 Free data ?/'ij[Re(A(2))/r,Im(A(i))] 

The discussion of the asymptotic structure of solutions of equation (|12p if the free data is given 
by ■i/'y [Re(A'^^^)/r, Im(A*^^-')] is not directly covered by the techniques of ^51. The reason for this is 
the appearance of terms with even powers of l/r in the source terms D^^ij. As it will be shown, 
these cases — which are the ones of more relevance for the present work — can be analysed under 
further assumptions on the free data. 
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The analysis in section[5]shows that free data of the form tfjij[A, J, QJ+^/jy [Re(A'-^-'),Ini(A''^-')/r] 
does not contribute to the regularity condition discussed in the present article. In order to 
have a non-trivial contribution from the second fundamental form in the regularity conditions at 
the critical points X^ one has to consider free data of the form ■!/'ij[Re(A'-^-')/r, Im(A^^-')]. The 
complications concerning this case will be discussed in the sequel. 

In a neighbourhood Ba of infinity, the form that the conformal metric acquires in normal 
coordinates can be used to decompose the differential operator appearing in equation p2|) , namely 

1 
3^ 



Lhv' = DkD^v' + ^D'Dkv" + r\v'' 



in the form 
where 



Ln + L 



+ -Lih, 



LoUi = d^dkv^ + -did^Vj, 



with A-'j and B •'^ functions of the metric coefficients and their first and second derivatives. They 
are smooth functions and satisfy 

Af, =0(r), B'''^ = Oir^). 

Properties of the operators Aq and L^. As before, let Vm denote the real linear space of 
homogeneous polynomials of degree m. Likewise, let Tim denote the space harmonic polynomials 
of degree m. If atj^.-.i^x^^ ■ ■ ■ cc'" £ Pm, cti-^.-.t^ = a(i^...i^), then it is an harmonic polynomial if 
and only if ai-^...i^ is trace free. It is well known that the space Vm can be written as a direct 
sum 

I m ^ "m ffi f "m— 2 ffi f ii-rn—A ® ' ' ' ■ 

Let s e Z, 

Ao = a'^Sfc, Ao : r^n. -> r^'^Vm, 

defines a bijective linear map if either s > or s < 0, \s\ is odd and to + s > — see [5]. If one 
wants to discuss the bijectivity of the Laplacian for functions of the form r'^Pm with s < 0, \s\ 
even and pm G Vm, one has to restrict the domain and the range sets. To this end define 

Elaborating from the proof of lemma 3 in 5 one finds the following lemma. 
Lemma 8. Let s G Z, with s < 0, \s\ even, then 



^^0 • ' ' 771, S ^ ' ' '> 



m,s 



is a bijective mapping if m + s > 0. 

The operator Lq has nice properties with regard to the spaces Qm. Indeed, if s G Z, then 

is a bijective linear mapping again if s > or s < 0, \s\ is odd and m + s > 0. As in the case of 
the Laplacian, in order to obtain bijectivity for s < 0, \s\ even, one has to restrict both domain 
and range. Let 

Q,n,s = {v' G C°°(M^R3) I v' G rm,s, v'x, = r\, V G P„,,4 , 

then again, following closely the arguments of [5] one can prove the following lemma — cfr. lemma 

11 in [g. 

Lemma 9. Let s G Z, with s < 0, \s\ even, then 

Co:r'Qm^s^r'-^Qm,s 
is a bijective mapping if m + s > 0. 
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Concluding the analysis. Returning now to the behaviour of solutions to equation (|12p with 
seeds V'ij[Re(A^^-')/r, Im(A''^-')], lemma [7] renders that 

£)*$y-[Re(A(2))/r,Im(A(i))] = r-^Sj[Rc{X'-^^)/r,lm{X'-^^)], S'j[Re(A(2))/r,Im(A(i))] = 0{r^), 

with S'j[Re(A(2))/r,Im(A(^))] G C°°{S). For the ease of notation, in what follows the affix 
[Re(A(^^)/r, Im(A*^^-')] will be dropped. One can write for some arbitrary m e N 

m 
k=3 

with s^ — o{r"^). The latter suggests considering 

m 
k—s 

SO that 

\/c=3 / \fe=3 / 

Using lemma 12 in reference [5] one readily finds that 

(ra \ jn m 

H ^'(fe) ) ^^2^1^ Hk) ) = '^"' I] "(fc) ' 
fe=3 / fe=3 k=3 

with Ujj,x — 0{r'^^^), U(^j G Qfc(;Ba). The functions Mk^ can also be written as 



"(fc) 



j=fc+3 



SO that 



m ?n 



^mE^w ^^ 'E E "k.)+^ 'E^wT^' 

\k^3 / k^3 j^k^S k^3 

7n j — 3 m 



'• 'EE"K.)+^ 'E" 



-2 

(fc)-R- 
j=6 fe=3 fc=3 



The above analysis suggests calculating the polynomials w^ s recursively if one can find solu- 
tions to the equations 



Lo(«w)=^"Mfe)' 6<fc<8 (62) 

I^o(^w)-^"MsW-E"Ufc))' ^<k<m. (63) 



3=3 

From lemma [9] one sees that one can only find polynomials i'/j,> solving ([62]) and ((63|) if the right- 
hand sides are in r^^Q^^s. Thus, in general, solutions to equation (fT2| with seed V'ij[Re(A(^')/r,Im(A'-^^)] 
cannot be expected to be smooth at i, and its asymptotic expansions may have terms involving 
the function Inr. 
The conditions 

slk) e Qk^3, Q<k<8, 

fe-3 



Sn. 



(k) 



E"0-.fc) ^2m, 9<fc<m, 



j=6 
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can, in principle, be reformulated as conditions on the free data ft.y and Re(A'-^'')/r, Im(A'^^), so 
that they cannot be independent of each other. We shall neither be concerned with a detailed 
analysis of these conditions nor with an explicit formulation of them as this goes beyond the 
scope of the present work. Through out it will be assumed that they are satisfied. Nevertheless, 
it is important to remark that there exists a class of free-data, (hij^ipij), with hij smooth for 
which this is the case: axial symmetric data — see |4]. The proof that this is the case is done 
by methods different to the ones discussed here and exploits that in axial symmetry, one has an 
explicit solution, ipij , to the momentum constraint given in terms of a potential. 

Once the polynomials v't^s have been determined, one is left to analyse the behaviour of the 
remainder v\). The remainder satisfies the equation 




Using the methods of theorem 17 in [5] one concludes that the right-hand side of this equation is 
in C""'^*~^'"(Sa). Standard results on elliptic regularity then yield w|j G C™'^'*'"(;Sa). Since m is 
arbitrary, one can conclude, by using a further argument that w]^ G C°°{Ba). 

Summarising, if w'[Re(A*^^-')/r, Im(A'^^^)] is the solution to equation (fT2| with with seed given 
by V'»j[Re(A(2))/r,Im(A(i))], then 

z;'[Re(A(2))/r,Im(A(i))] = i;i[Re(A(2))/r,Im(A(i))] + i;^[Re(A(2))/r,Im(A(i))], 
with 
^i[Re(A(2))/r,Im(A(i))] e Qoo(S„), t;i[Re(A(2))/r,Ini(A(i))] = 0{r^), v\[^e{\^^^)/r,l^{\^^^ e C^{Ba) 

B.2.7 The condition V'^iV'*-' G E°^{Ba) 

In order to be able to use theorem [7] with the ijjij which has been constructed in the previous 
paragraphs, one has to verify that 

V'., = ^y[A, J,Q] +7A,,[A(i)] -f V^,,[A(2)] + {Cv),,[A,J,Q] + (£„«),j[A(i)] + {Chv\,[\^% 

satisfies r^ipijip^^ G E°°{Ba). A lengthy calculation using the ideas of section 4.3 in [S] together 
with lemma [7] shows that this is the case. 
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